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A COMPUTATIONAL SYSTEM FOR 
AERODYNAMIC DESIGN AND ANALYSIS OF 
SUPERSONIC AIRCRAFT 

PART 3 - COMPUTER PROGRAM DESCRIPTION 


W. D. Middleton, J. L. Lundry, and R, G. Coleman 
Boeing Commercial Airplane Company 


1 . 0 SUMMARY 


An integrated system of computer programs has been developed for 
the design and analysis of supersonic configurations. 

The system consists of an executive driver and seven basic 
computer programs including a plot module, which are used to build 
up the theoretical force coefficients of a selected configuration. 

Documentation of the system has been broken into three parts: 

Part 1 - General Description & Theoretical Development 

Part 2 - User’s Manual 

Part 3 - Computer Program Description 

This part, the computer program description, contains schematics 
of the program structure and written descriptions of the 
individual overlays and subroutines. 

These three documents supersede NASA contractor reports CR-2520, 
CR-2521, and CR-2522 which described an earlier version of the 
system. 

Interactive graphics for use with the system are optional, 
employing the NASA-LRC CRT display and associated software. 

The computer program is written in FORTRAN XV for a SCOPE or 
KRONOS operating system. It is designed for the CDC 6600 series 
of computers and is overlay structured. The system requires 
approximately 115OOO 0 (octal) central memory words and uses eight 
disc files in addition to the input and output files. 




2.0 DISCUSSION 


A schematic of the design and analysis system overlay structure is 
shown in figure 2.0-1. The system is a single overlaid program, 
with the executive driver as the main overlay and the basic 
programs as primary overlays. The basic programs manipulate input 
(geometry module) , draw a picture of the configuration (plot 
module) , or perform design or analysis calculations. 

The format of the computer program documentation is to present 
schematics or block diagrams of the major program structure, 
together with subroutine descriptions, for each module developed 
under the design and analysis system contract. The plot and far- 
field wave drag modules are not included in this procedure, since 
they are described in other NASA documentation (references 1 and 
2 ) - 

The description of the overlay structure follows the convention of 
labeling overlays with octal numbers, but calling them with their 
decimal equivalents. 

A typical test case and associated program output is given in the 
User’s Manual (part 2) . 


File Usage 


File usage in the system is assigned as follows: 


1 

2 

3 

5 

6 

9, 10, 11, 12 
9, 11, 12 


basic geometry storage 
interface file 
restart data 
card input 
output 

storage files for far- field wave 
drag module 

storage files for plot module 


Program Structure 

A block diagram of the design and analysis system is shown in 
figure 2.0-2. The largest element of the system with the NASA-LRC 
graphics software attached occurs with the geometry display module 
loaded, and is approximately 115000 B . A "stripped” version of the 
system (graphics software excluded) would have its largest core 
requirement with the wing design module loaded, and would be 
approximately 77000 8 . 

These core sizes are for an absolute version of the program. 
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_ __ _ _ _ _ ... awaa , a . 1 Drag-due-to-lift 

Zero Lift Wave Drag (Wave and vortex) 


FIGURE 2.0- 1. -INTEGRA TED SUPERSONIC DESIGN AND ANALYSIS SYSTEM 









FIGURE 2.0-2 -SYSTEM 0 VERLA Y STRUCTURE 




2.1 EXECUTIVE MODULE 


The executive level ( 0,0 overlay) is used to read executive 
control cards and request execution of -tdie basic programs (primary 
overlays) as instructed. The executive cards are described in the 
user’s manual (part 2) and summarized under subroutine CHECKIN 
(page 6 ) . 

Program A3 8 9 AO 0 

PURPOSE: To read program executive cards and call primary 

overlays. 

METHOD: A389A00 is the program name assigned to the 

executive level (A389, 0, 0). It reads executive 
cards, calls subroutine CHECKIN to find the 

corresponding overlay number, calls the geometry 
module to read or sort input, then calls the 

appropriate primary overlay for problem execution. 

CRT is a special executive card used only to turn 
on (or off) the graphics routines as described in 
the user's manual. 

INPUT: Executive cards (see user's manual) 

SUBROUTINES 

CALLED: CHECKIN 


Program CHECKIN 

PURPOSE: To identify primary overlay number corresponding to 

executive card. 

METHOD: The overlay number for the different primary 

overlays is determined by variable PROG, in common 
block/SAVl /. CHECKIN sets the value of IPRG by 
finding the executive card word corresponding to 
PROG. The correspondence is: 


EXECUTI VE_CARD (PROG) IPRG 


PL0T 2 

FFWD 3 

SKFR 4 

GE0M 5 

NFWD 6 

ANLB 7 

WDE8 8 

FSUP 9 

WGUP 10 
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USE 


CALL CHECKIN 


INPUT: Variable PROG in common block SAV1 

SUBROUTINES 
CALLED: None 


Function TBLU1 
Integer Function L00KUP 

PURPOSE: To perform linear or second order interpola- 

tion from one-dimensional array, 

METHOD: TBLU1 and the associated integer function L00KUP 

are general purpose interpolation routines, for 
either linear or second order interpolation. The 
program call is; 

% = TBLU1 (X, XX, Y, MD, N) 

where 

X independent variable 

XX array containing independent 

variable 

Y array containing dependent 

variable 

MD code defining interpolation type 

1 - linear 

2 = second order 

N number of variables in array 

XX or Y 

B dependent variable (answer) 

Array XX must be monotonically increasing. 

USE: B =TBLU1 (X, XX, Y, MD, N) 

INPUT: As described above 

SUBROUTINES 

CALLED: TBLU1 calls associated subroutine L00KUP 
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2.2 GEOMETRY MODULE 


The geometry module is primary overlay 5. It contains subprograms 
to read, store, or display the configuration geometry, and to set 
up the input for the other basic programs. A schematic of the 
geometry module is shown in Figure 2.2-1. 


Program PDPACK 

PURPOSE: To route geometry handling requests from the 

executive to geometry routines, based on executive 
control cards. 

METHOD: Program PDPACK is the primary level of the geometry 

module. It is entered from the executive to store, 
add, or change input data; to update wing camber 
surface or fuselage basic geometry; or to enter the 
interface routines to set up input data for the 
other modules. A schematic of PDPACK is shown in 
figure 2.2-2. 


The executive control cards corresponding to 
variables in PDPACK are as follows: 


PRINCIPAL 


EXECUTIVE 

PDPACK 

PROGRAMS OR ' 

CARD 

VARIABLE 

SUBROUTINES CALLED 

GE0M NEW 

IFIRST=0 

EDITS 


GE0M 

IFIRST=1 

INPTS, 

EDITS 

PL0T 

IPRG=2 

INPTS, 

GE0MPLT 

FEWD 

IPRG=3 

INPTS, 

GE0M8O 

SKFR 

IPRG=4 

INPTS, 

GE0M158 

NEWD 

IPRG=6 

INPTS, 

GE0M916 

ANL2 

IPRG=7 

INPTS, 

GE0M2O 1 

WDE2 

IPRG=8 

INPTS, 

GE0M253 

FSUP 

IPRG=9 

INPTS, 

FUSUPD 

WGUP 

IPRG=10 

INPTS, 

NEWCAM 


Configuration geometry is read from cards in 
program EDITS and stored on tape 1 when the 
geometry module is not in core. Program INPTS 
retrieves the geometry from tape 1 when the 
geometry module is called. Subroutines WRGE0M and 
PRINTG write the configuration geometry on tape 1 
for storage or print the geometry onto the output 
tape, respectively. 


8 



f~ — — 

i Executive 


1 


GE0M or 
GE0M NEW card 


I 

L 


Read control card 

— ~i 


I 

._l 

SKFR, ANLZ, etc. 


Read 

or 

update 

basic 

configuration 

geometry 


Store basic 
geometry 
when not 
in core 


Interface subroutines 
(Set up input for selected program) 




Near- 

Far- 

Wing 


Plot 

Skin 

field 

field 

design 

Lift 

friction 

wave 

wave 

and 

analysis 



drag 

[_____] 

drag 

optimization 





Read additional input as necessary 

i 1 

Write input tape for selected program 



(Basic geometry + additional input, 
structured as required) 


1 


Basic 

geometry 

file 

Uape 1) 


t 


Input file 
for selected 
program 
(tape 2) 


FIGURE 2.2-1. -SCHEMATIC OF GEOMETRY MODULE 





Exit 


Exit 




CRT 

display 


FIGURE 2.2-2.-PR0GRAM PDPACK SCHEMA TIC 




USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


The interface programs write data onto tape 2, 
which becomes the input tape for the individual 
basic modules. 

CALL OVERLAY (GEOM # 5, 0, 0) 

Executive card variables listed above. 


See schematic on page 10 . 


Subroutine WRGE0M 

To write geometry data onto tape 1 for storage , or 
onto tape 2 for plot program or far- field wave drag 
program. 

Subroutine WRGE0M uses the input format of the 
NASA-LRC plot program and is used to write the 
basic geometry data onto tape 1 when the geometry 
module leaves core, or to write the geometry data 
onto the interface tape (2) as a part of the plot 
or far-field wave drag program interfaces. 

CALL WRGE0M 

Configuration geometry (see user’s manual) 


None 


Subroutine ZF0D 

To calculate the Z distance between the nacelle 
centerline and the local wing camber surface. 

The basic configuration geometry input (read in 
EDITS) allows the nacelles to he input with either 
of two Z dimensions: in the Z system of the 
configuration coordinates, or with Z measured to 
the local wing camberline. Subroutine ZP0D is 
used, regardless of which 8 definition is input, to 
calculate the other definition and store it in the 
basic geometry. 

If the nacelle origin is forward or aft of the wing 
planform, the Z distance is measured from the wing 
leading edge or trailing edge (whichever is closer) 
at the same Y station. 
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USE: 

INPUTS: 


SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


As a special case, 2P0D is also used to compute the 
distance between the fuselage centerline and the 
wing camber line at the side of the fuselage (in 
program GE0M916) . 

CALL 2P0D (NN) 

NN Nacelle origin index 

P0D0RG Nacelle origin data 

XAF, 

WAF0FG , Wing definition (plan form, 

WAF0RD, thickness, camber line) 

TERD 


None 


Subroutine PRINTG 

To write conf iguration geometry onto output tape 

( 6 ) . 

After conf iguration geometry is read or changed, 
subroutine PRINTG is used to write the current 
geometry definition* In the case where the 
fuselage or wing camber surface is updated 
(executive cards FSUP or WGUP) , only the fuselage 
or camber surface is output. This is controlled by 
variable ISKIP1 (=98 for fuselage only, =99 for 
camber surface only) . 

CALL PRINTG 

Geometry definition (see user’s manual) 

I SKIP 1 


None 


Subroutine NEWCAM 

To interpolate a wing camber surface for basic 
geometry definition (T20RD) from wing design or 
analysis definition (WE0RD) 

The camber surface definition generated or used by 
the wing design or analysis modules, called W80RD, 
may be different from the basic geometry definition 
(T20RD) . This can happen either because the 
definition is being created by the design program. 
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USE: 

INPUT: 


or because the analysis program input option for W2 
0RD was used* NEWCAM linearly interpolates for T2 
0RD from W20RD, and is called by program PDPACK if 
executive card WGUP is read. NEWCAM is also 
automatically called by the analysis program 
interface if the W20RD option is used. 

CALL NEWCAM (FACT0R) 

FACT0R Multiplier used to scale camber 
surface ordinates, if desired. 

N0PCT, 

P, 

J BYMAX # Camber surface definition WZ0RD 

Y, contained in common block /CAMBER/ 

W20RD, 

IF2C 


XAF, Camber surface definition (T20RD) 

WAF0RG, of basic geometry contained in 

T20RD common block /WING/ 


SUBROUTINES 

CALLED: INTERP 


PURPOSE: 

METHOD: 


USE: 


Subroutine INTERP 

To interpolate two-dimensional array. 

INTERP is used to linearly interpolate W20RD 
definition (used in conjunction with subroutine 
NEWCAM) . 

CALL INTERP (A, 0, N, X, Y, M) 


INPUT 


A First array location for abscissa 

values 

0 First array location for ordinate 

values 

N Number of values in A and 0 arrays 

X Array of abscissa for ordinate 

outputs 

Y Array of calculated ordinates 

M Number of values in X and Y arrays 

output 


SUBROUTINES 
CALLED: None 
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PURPOSE: 

METHOD:. 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


Program GE0MPLT 

To write input tape for PLOT module. 

GE0MPLT is the interface program between the 
executive and the PLOT module, entered when the 
executive card PL0T is read. The interface reads 
the plot title card and configuration codes, then 
uses WRGE0M to write the corresponding basic 
geometry data onto tape 2 (the program input tape) . 
GE0MPLT then reads the plot view cards and writes 
them onto tape 2. 

CALL 0VERLAY (GEOM, 5, 2, 0) 

See user's manual for input description 


WRGE0M 


Program GE0M8O 

To write input tape for far-field wave drag module. 

GE0M8O is the interface program between the 
executive and the far-field wave drag program, 
entered when the executive card FFWD is read. The 
interface reads a title card and the configuration 
codes, then uses WRGE0M to write the corresponding 
basic geometry data onto tape 2. GE0M8O then reads 
the case and restraint cards (if used) and writes 
them onto tape 2. 

CALL 0VERLAY (GEOM, 5, 3, 0) 

See user's manual for input description 


WRGE0M 


Program GE0M1 58 

To write input tape for skin friction module. 

GE0M158 is the interface pregram between the 
executive and the skin friction drag module, 
entered when the executive card SKFR is read. It 
reads title information, flight conditions, and 
configuration data, and structures those inputs 
together with basic geometry to write the input 
tape (2) for the skin friction module. 
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USE: 


INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


CALL 0VERLAY (GEOM, 5, 4, 0) 

See user’s manual for input description 


None 


Program GE0M916 

To write input tape for near-field wave drag 
module. 

GE0M916 is the interface program between the 
executive and the near-field wave drag module, 
entered when the executive card NFWD is read. It 
reads title information, Mach number, etc. , and 
structures those inputs together with basic 
geometry to write the input tape (2) for the near- 
field wave drag module. 

The near-field program considers the conf iguration 
to be uncambered; the wing and fuselage have 
thickness but have a flat mean-line. Wing height 
in the side of the fuselage is preserved, however, 
being the distance from the fuselage centerline to 
the wing at the side of the fuselage. In GE0M916, 
the wing height dimensions are computed by a 
special use of 2P0D, in which the fuselage 
centerline is treated as a nacelle origin in a 
series of calculations. 

The nacelle 2 dimension used in the near-field 
program is the distance from the local wing 
camber line. 

CALL OVERLAY (GEOM, 5, 6, 0) 

See user’s manual for input description. 


2P0C, NEWAREA 


Subroutine NEWAREA 

To check for and remove "steps" in input fuselage 
area definition corresponding to inlet or exit 
stream tubes. 

The basic geometry fuselage definition allows for 
input of the fuselage in four segments. If the 
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first: cross-section of segment 2 is input at the 
same X station as the last cross-section of segment 
1, and the cross-sectional areas of the two inputs 
are different, then the far-field wave drag program 
extends a streamtube to account for the area 
"step” , so that no area discontinuity occurs in the 
wave drag calculations. This can occur, also, 
between other segments. The near-field wave drag 
fuselage calculations cannot accommodate an area 
step at a fuselage station other than the first or 
last. If one occurs, subroutine NEWAREA removes 
the step by collapsing the fuselage to a solid body 
having the same area growth. 


USE: 

INPUT: 


SUBROUTINES 

CALLED: 


NEWAREA also reduces the fuselage area definition 
to 50 values in X, 8 and area if more than 50 total 
values were input in all fuselage segments, since 
the near-field program allows only 50 points. 

CALL NEWAREA (X, 8, A, L) 

X Output array of fuselage X values 

for analysis program. 

8 Output array of fuselage camberline 

8 values. 

A Output array of fuselage cross- 

sectional areas, 

L Number of values in X, 2, or A arrays 

/B0DY/ Common block containing fuselage 

basic geometry definition- 


Ncne 


Program GE0M2O1 

PURPOSE: To write input tape for lift analysis module. 

METHOD: GE0M2O1 is the interface pregram between the 

executive and the lift analysis module, entered 
when the executive card ANL2 is read. It reads 
title information, Mach number, etc. , and 
structures those inputs together with basic 
geometry to write the input tape (2) for the lift 
analysis module. 


The lift analysis interface allows the wing camber 
surface to be read as the W80RD definition, or 
passed in through common block /CAMBER/, in 


16 



USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 


addition to the basic geometry definition. If W2 
0RD is used, the basic geometry is automatically 
updated by means of subroutines NEWCAM and WRGE0M. 

In addition, since the lift analysis program logic 
does not permit step discontinuities in fuselage 
area between the most forward and aft X stations, 
subroutine NEWAREA is used to prepare the fuselage 
cross-sectional area definition. (NEWAREA is 
described in connection with GE0M916). 

The nacelle origin 2 dimension is the vertical 
distance to the local wing camber line. 

CALL 0VERLAY (GEOM, 5, 7, 0) 

See user's manual for input description 

NEWAREA, NEWCAM, WRGE0M 


Program GE0M253 

To write input tape for wing design module. 

GE0M253 is the interface between the executive and 
the wing design module, entered when the executive 
card WDE2 is read. It reads title information, 
Mach number, loadings data, etc. , and structures 
those inputs together with basic geometry to write 
the input tape (2) for the wing design module. 

The wing design program may employ pressure fields 
generated by the near- field wave drag module or 
analysis module. The existence of these pressure 
fields is tested in the corresponding common 
blocks. This involves the cedes B0DCPX, B0DUPX, 
and CPNACX . The actual use of these pressure 
fields is controlled by input of the corresponding 
loading numbers. 

The nacelle 8 origin is the vertical dimension to 
the local wing camber surface. 


The restart data deck is substantial. Logic for 
this part of the interface is copied from the wing 
design program input. 

CALL 0VERLAY (GEOM, 5, 8, 0) 

See user's manual for input description 
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SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 


RDWRT, RDWRTE 


Subroutine RDWRT, RDWRTE 

To read input data cards and write the data on the 
interface file (2) in the format specified. 

If subroutine RDWRT is called, data is read using 
FORMAT (10F7. 0) and written using FORMAT { 1 0F1 0. 4) , 
If subroutine RDWRTE is called, data is read and 
written using FORMAT (4E20- 1 3) • 

CALL RDWRT (NUM) 

CALL RDWRTE (NUM) 

NUM Number of data values to transmit. 


None 


Program FUSUPD 

To update the fuselage definition in the basic 
geometry to the optimum area distribution generated 
by the far -field wave drag program. 

The far-field wave drag program contains a fuselage 
optimization feature, which area-rules the fuselage 
subject to input constraints. FUSUPD uses the 
optimized definition to update the basic fuselage 
geometry, and is accessed by means of the executive 
card FSUP, 

FSUP first reads the value of 0PHOW, which controls 
the optimization interpolation, (0PH0W = -, 
interpolate at original fuselage X stations; 0PH0W= 
+ , interpolate at 50 equally spaced X stations) . 
It then changes the basic geometry definition to 
the optimum body area distribution, contained in 
common block/0PBOD/, The original fuselage 2 
definition is preserved. Fuselage perimeters are 
proportioned to the new area distribution. 


CALL 

0VERLAY (GEOM, 5, 9, 

0) 

OPH0W 

Interpolation code 

X0P 

Optimum fuselage 

definition- X0P 

A0P 

array, A0P = area 

array, J0P = 

J0P 

number of X0P or 

A0P values. 
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SUBROUTINES 

CALLED” 


PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 


PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 


TBLU1 


Program INPTS 

To read basic geometry data frcir tape into core. 

Program INPTS is used to reread the configuration ' 
basic geometry back into core when the geometry 
module is accessed. INPTS uses the same format as 
program WRGE0M, which was used to store the basic 
geometry on tape 1. 

CALL 0VERLAY (GEOM, 5, 10, 0) 

Tape 1 


None 


Program EDITS 

To read configuration geometry from cards. 

Program EDITS is used to read the configuration 
geometry from cards, and is accessed by either the 
GE0M NEW or GE0M executive cards. {GE0M NEW zeroes 
all the configuration J1, J2, etc., codes, so that 
EDITS reads all new data; GE0M preserves all 
existing codes, so that new geometry read replaces 
or adds to existing geometry read by INPTS) . 


The input format of program EDITS is basically that 
if the NASA-LRC plot program, with a few additional 
input variables added as noted in the user’s 
manual . 

CALL 0VERLAY (GEOM, 5, 11, 0) 

See user's manual for input description 

SP0D 
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2.3 SKIN FRICTION MODULE 

The skin friction module is primary overlay 4, Its principal 
subprograms calculate the wetted areas of the configuration, and 
then the skin friction drag for a series of input flight 
conditions. A schematic of the skin friction module is shown in 
figure 2.3-1. 

Program TEA-158A 

PURPOSE: To read configuration data, flight conditions, and 

calculate wetted areas and lengths to be used in 
skin friction drag calculations . 

METHOD: Program first reads all input data, which consists 

of configuration geometry and either Mach number- 
altitude or Mach number-Reynolds number flight 
conditions. 

For each configuration component, the corresponding 
wetted areas and reference lengths are computed. 
In the case of the wing or canard (s) or f in (s) , the 
parts are broken into strips to permit more 
accuracy in determining an average skin friction 
coefficient. (The wing is broken into 
approximately 50 strips, canard or fin into 10 
strips.) 

The wetted areas and references lengths are passed 
to subroutine DRAG for the skin friction 
coefficient calculations. 

USE: CALL OVERLAY (SKFR, 4, 0) 

INPUT: Input is read from the interface tape (2) , set up 

by the skin friction interface program in the 
geometry module. 

SUBROUTINES 
CALLED: DRAG 


Subroutine DRAG 

PURPOSE: To calculate configuration skin friction 

coefficients and print answers for all input flight 
conditions . 

METHOD: Given the set of configuration wetted areas and 

input flight conditions, the skin friction 
coefficients for each component are computed and 
summed to obtain the total skin friction drag. If 
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• Reads input 


• Displays 
results 
on CRT 



• Computes wetted areas 
and reference lengths 


• Calculates friction drag 
for input flight conditions 


• Skin friction coefficient 
from T' method 


• Standard atmosphere data 


FIGURE 2.3-1. -SKIN FRICTION MODULE SCHEMATIC 



the graphics routines are activated, the friction 
drag coefficients are also printed on the CRT. 

USE: CALL DRAG 

INPUT: Input is passed in common blocks from program TEA- 

1 58A 

PRINCIPAL 

VARIABLES: Same as program TEA-158A, plus: 

SCAM Scale factor to convert configura- 

tion reference lengths to feet for 
use in subroutine FICT 
SREF Wing reference area 

SWETT Total configuration wetted area 

CDFT Total skin friction drag coefficient 

AM, Mach number - altitude flight 

AL, condition input AM = Mach number, 

DELT, AL = altitude (feet) , DELT = 

SCAM0D temperature deviation from standard 

(°F) , SCAM0D = input scale factor. 

AM, Mach number - Reynolds number 

RNPFL, condition input. RNPFL = Reynolds 

SCAM0D, number per foot/1,000,000 T0TEM = 

T0TEM total temperature (°Rankine) 

SWETRB Fuselage wetted area 

FUSL Fuselage reference length 

SWETB Fuselage wetted area corrected for 

overlap areas of wing, fin(s), and 
canard (s) 

AREAJ Planform area of a wing strip 

WINGL Length of wing strip 

LCOUNT Total number of wing strips 

SWETW (nominally 50) total wing wetted 

area 

TSWTNA Wetted area of nacelle (s) at first 
TP0DL input origin, corresponding length, 

SWETNA total wetted area for all nacelles 

NP0D Number of nacelle origins 

SWTFN Planform area of a fin strip (10 

CHDFN total strips) , corresponding strip 

NFIN length, number of fins 

TSWTC Planform area of canard strip 

TCCAN Corresponding strip length 

NCAN Number of canards 


22 



SWETXP Wetted area of an arbitrary con- 

RLXP figuration part, corresponding 

NXTPT reference length, number of 

arbitrary parts, 

OUTPUT: Output consists of flight conditions, wetted 

area and skin friction coefficient buildups, 
total skin friction drag coefficient and con 
figuration wetted area. 

SUBROUTINES 

CALLED: FICT, DISP158 


Subroutine FICT 

PURPOSE: To calculate the turbulent skin friction 

drag coefficient for a given reference length, 
Mach number, Reynolds number, and total 
temperature. 

METHOD: The skin friction coefficient is calculated 

from the T* method described in the aerody- 
namic theory document (part 1) . 

USE: CALL FICT (AMX, ALX , EEL, C) 

INPUT: AMX Mach number 

ALX Altitude, feet 

EEL Reference length 

PRINCIPAL 

VARIABLES: TI Total temperature, degrees Rankine 

SCAM Scale factor to convert input 

reference length to feet 
RI Free stream Reynolds number 

OUTPUT: C Skin friction coefficient 

SUBROUTINES 

CALLED: ATM062, SKIN 


Subroutine ATM062 

PURPOSE: To provide standard atmospheric data 

METHOD: Program uses 1962 standard altitude definition 

(reference 3) . 

USE: CALL ATM062 (Z* TEMP, SIGMA, AX) 
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INPUT: 

OUTPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 

OUTPUT: 

SUBROUTINES 

CALLED: 

ERROR 

RETURN: 


Z Geometric altitude, feet 

TEMP Temperature, degrees centigrade 

SIGMA Density ratio 

AX Speed of sound, knots 

None 

Subroutine SKIN 

To iterate for skin friction coefficient. 

Program is used to solve the Karman- Schoenherr 
equation 


using 


« 242 
VCF 


- Ion 


10 


(CF * KEY) 


value for CF. Solution is iterative, and 
is satisfied when successive iterations agree 
within .0001 per cent. A maximum of 50 
iterations is allowed. 


CALL SKIN (REY, CF) 

REY Reynolds number 

CF Skin friction coefficient 


None 


Program uses 50th iteration for CF if 
convergence does not occur, and prints error 
message. 
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2.4 NEAR-FIELD WAVE DRAG MODULE 




The near- field wave drag module is primary overlay 6. It contains 
principal subprograms to calculate the near- field pressure data 
and drag coefficients and to display the calculated results. 
Schematics of the principal program structure are given in figures 
2-4-1, 2.4-2, and 2.4-3. 

Program TEA-356 

PURPOSE: Near-field wave drag primary overlay 

METHOD: Program TEA-356 is the primary level of the 

near-field wave drag program. It calls the 
input overlay and contains the Mach number 
loop which calls the calculation overlays 
and the graphics display (overlay 3) . 

USE: CALL OVERLAY (NFWD, 6, 0, 0) 

INPUT: Executive card NFWD 

SUBROUTINES 

CALLED: Overlays 6,1 to 6,6 and primary level subroutines 


Subroutine NB0DY 

PURPOSE: To calculate the Whitham F (y) function for an 

arbitrary body of revolution, then compute surface 
pressure distribution, wave drag, and wetted area. 

METHOD: Subroutine NB0DY first computes the Whitham F (y) 

function, according to the Stieltjes integral 
formulation. The equations used are summarized in 
the theoretical document (part 1). The F (y) 
function is computed at TNCUT+ 1 X-stations (where 
TNCUT is an input in the main program, usually 50) . 
To remove minor irregularities in the final F (y) 
function, extra points are inserted between each of 
the X-stations, F(y) is computed there also, and a 
five point smoothing formula applied. 

The F{y) function is also computed at X-stations 
downstream of the last body station, to define the 
"tail” on the F (y) function which is needed in the 
area-balancing technique for pressure coefficient. 
The loop for the F(y) equations begins at statement 
60, which results in the unsmoothed F (y) array 
stored in TFY, with corresponding y values in BT0 , 
JFY total values. After smoothing and extending 
the F(y) function aft of the body, the final F (y) 
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function is stored in TFTAU, with corresponding y 
values in TTAU, JST0 total values. 

The body surface pressure calculations loop starts 
at D0 330, using the equations given in the theory 
document (part 1) . The body is again broken into 
TNCUT+1 segments, and the resulting thickness 
pressures stored in TCPBY, with corresponding X 
values in TXCPBY. The wetted area calculation and 
the integration of the thickness pressures to get 
wave drag is performed in loop £0 340. 

Subroutine NB0DY will handle bodies having open 
front or aft ends, but assumes smooth geometry in- 
between. 

USE: CALL NBODY (TXB, TPB, TNXB, TNCUT, TTAU, 

TFTAU, JST0, WAVE, SNWET) 

INFUT: TXB Body definition. TXB = x array, 

TRB, TRB = radii array, TNXB values in 

TNXB X or radius. 

TNCUT Number of body segments in F(y) 

and thickness pressure calculations 
TTAU F (y) function. TTAU = y, 

TFTAU, TFTAU = F (y) , JST0 values in y or 
JST0 F (y) . 

WAVE /C dA = wave drag/q 

SNWET Wetted area of body, assuming 

circular cross-sections 

SUBROUTINES 
CALLED: TELU1 


Subroutine NFWDR 

PURPOSE: To calculate influence coefficient function for 

wing thickness pressure calculations. 

METHOD: NFWDR calculates R function for a selected grid 

element according to equation given in theory 
document (Part 1) . 

USE: CALL NFWDR (NS1, DELTN1, DELTN2 , AN1) 

INPUT: NS1 X distance from field point element 

DELTN1 lateral distances froir field point 

DELTN2 element. By + .5 

AN 1 R function 
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None 


SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 


SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


Subroutine NINTPLT 

To generate a detailed definition of the tilted 
F(y) function for a selected right-running leg. 
Subroutine NPSIG identifies all right-running legs 
of the tilted F (y) function. NINTPLT is used to 
enrich the definition for use in later 
interpolations. The selected leg may either be a 
"base" leg, or a leg that possibly contains an 
area-balancing solution. The maximum number of 
points used to enrich the leg definition is 300. 
The enriched definitions end up in the following 
arrays: 


Tilted F (y) 

F(y) 

Integrated area 
under tilted F (y) 


Base leg 

BT0 

BFTAU 

BARE A 


Other leg 

CT0 

CFTAU 

CAREA 


CALL NINTPLT (SUM, TPRE, FPRE, NL, NR, NSET, 
TFTAU) 

SUM Area under tilted F(y) function 

corresponding to left end of leg 
TPRE, FPRE Value of tilted Y, (F (y) at left 
end of leg 

NL, NR Storage index of tilted F (y) array 

(arrays TT0 = tilted y, TFTAU = F (y) 
NSET Index for leg definition 

1 = base leg 

2 - succeeding leg 

TFTAU F (y) array 


None 


Subroutine NINTSC 

To locate the wing-fuselage intersection and 
compute corresponding exposed wing wetted area. 

The wing-fuselage intersection code is controlled 
by input ANYB0D. If positive, ANYB0D defines the 
inboard end of the wing for purposes of the 
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II 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


thickness definition, whether a fuselage is present 
or not. If ANYB0D is negative, the intersection of 
the wing and fuselage is calculated by solving for 
the point at which each constant per cent chord 
meanline of the wing crosses through the fuselage 
surface (assuming the fuselage to be circular) . 

The wing- fuselage intersection is stored in common 
block /WBINT/. 

WX X location of intersection 

WY,WE , VJT Corresponding Y,E, and thickness 
values 

CALL NINTSC (WSW) 

WSW Exposed wing wetted area 


TE.LU1 


Subroutine NPSIG 

To calculate the near- field pressure signature for 
a body of revolution at a given distance from the 
body. 

Given an F (y) function, Mach number, and radial 
distance from the body, NPSIG is used to calculate 
the resultant pressure signature. The method is 
described in the aerodynamic theory document (part 
1 ). 

The basic steps in the calculation sequence are as 
follows: 

1* The tilted F (y) function is constructed from 
the input F (y) function, and the running sum 
of area under the tilted F (y) function 
computed. (The input function is defined by 
arrays TTAU, TFTAU and the tilted function by 
TT0, TFTAU, with the area corresponding to TT0 
in TAREA) - 

2. The bow shock position is identified by the 
left-most value of TT0 for which TAREA - 0. 
(Interpolation of the TAREA array is used for 
accuracy. ) 

3. The tilted F(y) function remaining after the 
bow shock is located is searched to identify 
all right-running legs. 
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4. All right -running legs are compared for common 
TAREA values, which would identify possible 
shock wave locations- (Again, interpolation 
is used to enrich the TAREA definition, using 
subroutine NINTPLT.) 

5. Valid shock wave solutions are identified, 
based on the left-most criteria described in 
the theory document. 

6. The resulting tilted and area-balanced F(y) 
function is converted to pressure coefficient 
by means of the multiplier C0NST and stored in 
arrays BT0 and CT0 (x and pressure 
coefficient) . The BT0 and CT0 arrays are cut 
off if x becomes greater than input value YED. 

USE: CALL NPSIG (TTAU, TFTAU, JST0, R, YED) 


INPUT: 


SUBROUTINES 

CALLED: 


TTAU 

TFTAU 

JST0 

R 

YED 


F (y) definition, Y = TTAU array, 

F (y) - TFTAU array, JST0 values 

Radial distance from body centerline 
Largest x-value of interest for 
pressure signature 


NINTPLT, TBLU1 


Subroutine XFM 

PURPOSE: To convert input planforra geometry to program grid 

system* 

METHOD: The wing planform is represented in the near-field 

wave drag program as a set of recilinear elements, 
described in the theory document (part 1) . Given 
the number of semi -span rows (TN0N) used to define 
the right-hand wing, XFM interpolates the planform 
definition at the mid- point of each semi-span row 
for the leading edge and trailing edge x value. 
These x values are converted to program scale (by 
means of RATI0 = Ywing-tip * BETA/TN0N) and stored 
in arrays TXLE and TXTE, with XLE0 and XTE0 
defining the wing centerline grid points. 

Storage of wing pressure coefficients and surface 
slopes uses a space-conserving technique: The 

leading edge of a spanwise element row is stored 
immediately after the trailing edge of the adjacent 
inboard row. The index array for the row storage 
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USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


is J2, which totals the number of elements stored 
up to the leading edge of a selected row. The 
allowable total number of elements representing the 
right hand wing is 2750; if this total is exceeded 
for a given value of TN0N, the program reduces TN0N 
until the 2750 limit is not exceeded. 

CALL XFM 

Input geometry in common blocks 


TBLU1 


Subroutine PRECALC 

To precalculate a standard set of R values for 
thickness pressure calculations. 

PRECALC sets up an array of R values for grid 
elements up to 20 units forward or laterally from 
any field point element (tc save time in 

recalculating often-used values) . The resulting 
array is stored in common block /RSAVE/. 

CALL PRECALC 

None 


NFWDR 


Program P916 

To perform thickness pressure calculations in near- 
field wave drag module. 

Program P916 is the main program of the near- field 
wave drag module. It solves for the thickness 
pressure distribution on the surface of an 
arbitrary wing-fuselage-nacelle configuration, and 
integrates the pressures over the surface to obtain 
wave drag. 

The equations used in the thickness pressure 
solution are given in the theory document (part 1) . 
A consistent nomenclature is used, as possible, 
between the theory equations and the Fortran 
variable names. 


33 



The principal program logic is as follows: 

1) Calculate fuselage pressure field acting on 
wing. This is done for each selected semi- 
span station (TYB2) in loop D0 130. The 
resultant fuselage pressure field is stored in 
PBWG , which is printed in loop D0 150. 

2) Calculate wing thickness pressure field. This 

is done for the selected semi-span stations in 
loop D0 830. The semi -span element row is 
NS TAR and the associated x elements are LSTAR. 
For each NSTAR and LSTAR, the upstream region 
of integration contains elements located at N 
and LVAR. The local wing slopes are D2DX, 
contained in array T2C. The loop D0 450 
performs the influence function times slope 
calculation described in the theory document. 
The velocity potential (PHI) is then 

calculated, and differenced to get pressure 
coefficient (CP) . 

After all pressure coefficients are calculated 
at a given semi- span station, they are 

smoothed in loop 710, resulting in CPAVG. 
Calculation of the drag terms (pressure times 
wing slope) employs these variables names: 


DBWG 

fuselage pressure acting 


on wing 

DAVG 

wing thickness pressure 

DNAC 

nacelle pressure acting on 


wing. 


3) Interpolate wing thickness pressures for 
output. This is done at the selected semi- 
span stations in loop D0 770. The wing 
thickness pressures are then output in loop D0 
840. At the same time, the common block 
containing the wing thickness pressures, 
/PLIM/, is completed to make it self- 
sufficient for use in other modules, where: 

FYB2 per cent semi-span array, 

MYFR values 

TNMPC per cent chord array, NMPC 

values 

VAR (X, Y) thickness pressure array 

The wing thickness pressure plus fuselage 
pressure field is then summed and output in 
loop D0 860. 

Integrate drag terms spanwise to get total 
wing drag. This is done in loop D0 900. 
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USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 


5) Calculate wing thickness pressure on nacelle 
volume term (subroutine NTW0NAC*) 

6) Calculate fuselage pressure drag plus wing 
thickness pressure on fuselage drag terras. 
This is done in subroutine BDRAG and results 
in drag terms DBD and DW0B , respectively* 

7) Convert drag terms to coefficient form, output 
wing section drag coefficients, nacelle drag 
coefficients, and final drag summary* This is 
done starting at loop D0 920, using the 
following nomenclature: 

YB2 semi-span fraction 

CHORD local wing chord 

CAS wing section drag coefficient 

BEF section drag of fuselage on 

wing 

CNACS section drag of nacelle on 

wing 

CINT sum Of CAS, BEF, CNACS 

DFR fraction of total wing drag 

occuring at YB2 

CDAVG wing thickness drag coefficient 

CDB fuselage thickness drag coef- 

ficient 

CDBOW interference drag coefficient 

of fuselage on wing 

CDVJOB interference drag coefficient 

of nacelle on wing 
CDTOT total drag coefficient 

The nacelle drag coefficient summary is 
printed in subroutine ADDUM. Empennage drags 
are printed in subroutine ADEMP. 

CALL 0VERLAY (NFWD, 6, 1, 6HRECALL) 

MK Mach number index variable (from 

overlay 6, 0). Configuration data 
from overlay 6, 2. 


See schematic, figure 2.4-2* 


Subroutine ADDUM 

To sum nacelle drag terms, convert them to 
coefficient form (based on wing reference area) and 
write them onto the output file. 
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METHOD: 


Program NACL generates all nacelle thickness drag 
terms except for the wing on nacelle and nacelle on 
wing terms. These are saved in common block 
/NDRAG / and printed by ADDUM when the main program 
(P916) writes the final drag summary. The drag 
coefficient nomenclature is as follows: 

TNISOL isolated nacelle wave drag coef- 

ficient 

TFB0N interference drag coefficient of 

fuselage on nacelles 

TFN0B interference drag coefficient of 

nacelle on fuselage 
TNN0N interference drag of ether 

nacelles on nacelle 
TFIMAG “image" drag coefficient of 

nacelle on itself 

TF0IMG "image" drag coefficient of 

other nacelles 

TDW0N interference drag coefficient 

of wing on nacelles 

USE: CALL ADDUM (SUMM) 

INPUT: SUMM Sum of nacelle drag coefficients 

SUBROUTINES 
CALLED: None 


Subroutine AEEMP 

PURPOSE: To write and sum empennage drag terms. 

METHOD: Program EMPNN computes thickness pressure distribution 

and wave drag of canard (s) or fin(s), and stores drag 
coefficients in common block /EMPDRG/. Subroutine 
ADEMP is used to output the drag coefficients. 


USE: 

CALL 

ADEMP (SUME) 

INPUT: 

SUME 

Sum of empennage drag coefficients. 

SUBROUTINES 



CALLED: 

None 



Subroutine BDRAG 

PURPOSE: To integrate fuselage thickness pressures over 

fuselage to get wave drag, and calculate wing-on- 
fuselage interference drag term. 
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METHOD: 


USE: 

INPUT: 


SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


The fuselage is broken into TNCUT lengthwise 
strips, and 5 radial segments (per half section) 
per strip- The fuselage thickness pressures and 
wing pressures acting on the carry-over region 
covered by the fuselage are integrated over the 
fuselage segments. Wing pressures are transferred 
aft along Mach lines to obtain the point of 
application on the fuselage. 

CALL BDRAG (TNCUT) 


TNCUT 

TNMPC, 

TPWE, 

NMPC 

/PLIM/ 

TXCPBY , 
TCPBY, 

NCPB 


Number of fuselage strips. (Set 
at 50 in main program) . 

Wing pressure field extended along 
fuselage centerline. TPWE - pres- 
sure coefficient array, TNMPC = X 
array, NMPC = number cf pressure coef- 
ficient values 

Common block containing wing 
thickness pressures 
Fuselage thickness pressure array 
TXCPBY = X, TCPBY = pressure 
coefficient, 

NCPB values 


GRAB, TBLU1 


Subroutine DBCALC 

To calculate wing slopes for all grid elements 
representing wing planform. 

The wing thickness definition is contained in 
common block /THICK/, consisting of the array T20RD 
(thickness profile) , TPCT (per cent chord array, 
N0PCT values) . The Y stations corresponding to TH 
0RD are the planform input stations, in array YLED , 
D2CALC first interpolates the T20RD array linearly 
at the midpoint of each semi-span element, storeing 
the results in A2C. The ABC array is then 
interpolated chordwise to obtain the dz/dx slope 
for each element, storing the slopes in T2C. The 
chordwise interpolation may be either linear or 
quadratic depending upon an input code. (The T8C 
storage scheme is described under subroutine 
TRNSFM.) 

Inboard of the wing side-of-body Y value (input 
ANYB0D, if ANYB0D = +) , the wing slopes are set to 
zero. 
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USE: 


CALL DSCALC 


INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 


Wing planform, thickness, ANYB0D 
(described above) . 

TELU1 


Subroutine NGETNP 

To interpolate a nacelle pressure signature at a 
given wing semi-span station, for use in 
calculating nacelle-on-wing interference drag. 

The nacelle pressure field, TCPN, is a two- 
dimensional array (for per cent chord and span) . 
The wing program, P916, performs wing thickness 
pressure calculations at selected span stations. 
It is convenient to have a one -dimensional array, 
versus X only at that span station, for use in 
interpolating nacelle pressures for the nacelle-on- 
wing interference drag. 

NGETNP interpolates TCPN to obtain the required 
pressure signature, which is stored in TXNCW (per 
cent chord) and TCPNCW (pressure coefficient) . 

CALL NGETNP (YB2) 

YB2 Semi-span fraction 

/NPF/ nacelle pressure field definition 

common block. 


None 


Subroutine GRAB 

To interpolate pressure coefficient at a 
given planform location from the wing thick- 
ness pressure table or fuselage -on- wing 
interference pressure table. 

Linear interpolation for pressure coefficient 
froir a two-dimensional array. 

CALL GRAB (X, Y, XCP, JTELL) 

X X coordinate on wing planform 

Y Y coordinate on wing planform 
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SUBROUTINES 

USED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 


XCP 

JTELL 


VAR 

PBWG 


Interpolated pressure coefficient 
Variable to select interpolation 
array 

1 = wing thickness pressure 

2 = fuselage-on-wing pressure 
Wing thickness pressure array 
Fuselage- on- wing pressure array 


None 


Subroutine NW0NAC 

To calculate buoyancy effect of wing thickness 
pressures acting on nacelle. 

The wing-on-nacelle interference drag term, 
consisting of wing thickness pressures acting on 
the nacelle area distribution, is computed in 
NW0NAC. Wing thickness pressures are transferred 
aft along Mach lines to a series of nacelle first 
runs. The resulting interference drag terms, C_dA, 
are stored in array, TDW0N, with storage index 
corresponding to the nacelle origin number. 

CALL NW0NAC 

Wing thickness pressure definition (VAR) , configu- 
ration geometry passed in through common blocks. 


TBLUI 


Subroutine WGRID 

To define wing grid system data. 

WGRID is a utility routine, identifying wing grid 
element data at a selected spanwise element row. 

CALL WGRID (JBY, XLE, STE, LLE, LTE, ALE, ATE) 

JBY selected element row 

XLE wing leading edge coordinate 

XTE wing trailing edge coordinate 

LLE leading edge element 

LTE trailing edge element 

ALE leading edge element function (at LLE) 

ATE trailing edge element fraction (at LTE) 
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SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


None 


Program NACL 

To calculate pressure field of nacelles and 

associated drag terms, 

NACL is a major overlay of NFWD , which is used to 
calculate the nacelle thickness pressures, and all 
nacelle interference terms except the wing-on- 
nacelle and nacelle-on-wing terms. A schematic of 
the program is shown in Figure 2.4-3. The 

calculation sequence is as follows: 

1. Construct the pressure field of the nacelles 
acting on the wing. To do this, a series of 
wing semi-span stations is set up (in NYSET) . 
The program then calculates the composite 
pressure signature from all nacelles acting at 
those Y stations, for nacelles first below the 
wing, then for nacelles above the wing (if 
any) . Either glance or wrap nacelle pressure 
field option may be used. 

The nacelle pressure field is shrunk to 20 X 
and pressure coefficient values at each y 
station, stored in common tlock NPF: 

TXPN x array, per cent chord 

TCPN pressure coefficient array 

TYP y semi-span stations (per cent) , 

NYP values 

2. Calculate composite pressure signatures and 

inter f erence drag terms between nacelles or 
between nacelles and fuselage. The composite 
pressure signatures are applied to the area 
growth of the affected component to get the 
drag force. 

3. "Image” effects of the nacelles are 

calculated. If the nacelle is located next to 

the wing, the pressure signature from itself 
reflects off the wing and back onto the 

generating nacelle. Similarly, the reflected 
signature from other nacelles may cause an 
interference drag. Drag interference due to 
reflected pressure signatures are calculated 
separately from the direct effects under (2) . 
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use: 


CALL OVERLAY (NFWD, 6, 4,0) 


INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


Configuration geometry, passed to NACL in 
common blocks. 


See schematic, figure 2*4-3. 


Subroutine BPINT 

To integrate the nacelle (s) pressure signature over 
fuselage surface to obtain interference drag. 

The composite nacelle pressure signature (contained 
in XVAL, CPVAL, JNEXT values) is applied to the 
fuselage area distribution obtained in subroutine 

B SETUP • The process is repeated for each nacelle 
origin- The integral of Cp* area is called BF0RCE- 

CALL BPINT (ZNAC, BF0RCE, JNEXT) 

2NAC Nacelle vertical deminsion, rela- 

tive to wing 

BF0RCE Resulting pressure force 

JNEXT Number of values in nacelle pres- 

sure signature 


TBLU1 


Subroutine BSETUP 

To calculate fuselage area growth for use in 
integrating nacelle-on- fuselage interference drag 
term. 


The fuselage is broken into 3 pieces: forebody, 
mid-body region (wing part) , and aft body- The 
mid-body region covers the X interval of the wing- 
fuselage intersection. 


Each fuselage piece is broken into 50 segments, and 
the area growth associated with each segment stored 
in an array (for use later in subroutine BPINT) • 
The arrays are: 


FBX, FBDA 
BMX 
BMDAU 
BMDAX. 


fore body X and area 
mid-body X 

area growth above wing 
area growth below wing 
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use: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 


ABX, ABDA aft-body X and area 

CALL BSETUP 

Fuselage geometry contained in /FUSLG/ 
wing-fuselage intersection contained in 
/WBINT/ 

TBLU1 


Subroutine NC0MBIN 

To combine two nacelle pressure fields, for 
nacelles located above and below the wing, into a 
single composite pressure field- 

If nacelles are located both above and below the 
wing, subroutine NACL generates two nacelle 
pressure fields: TCPDN (nacelles below) and TCPN 

(nacelles above) . NC0MBIN sums these two nacelle 
pressure fields, to use in computing the thickness 
pressure interference drag- The resulting pressure 
field is stored in TXPN (per cent chord array) and 
TCPN (pressure coefficient array) , with NYP semi- 
span per cents contained in TYP 

CALL NC0MBIN 

Nacelle pressure fields described above 


TBLU1 


Subroutine NCUT0FF 

To delete portions of a nacelle pressure field. 

In constructing nacelle pressure signatures at an 
arbitrary distance from the generating body, 
subroutine NACL calls NCUT0FF to delete parts of 
the signature intercepted (or not reflected) by 
intermediate configuration components. The 

pressure signature is contained in BT0 AND CT® (x 
and pressure coefficients, respectively) . 

CALL NCUT0FF (XI, X2) 

Xl, X2 Bounding values of X, between which the 

pressure coefficients are set to zero. 
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SUBROUTINES 

CALLED: 


PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


TBLU1 


Subroutine NGLANC 

To perform cutoff of nacelle pressure signature 
according to glance feature of nacelle pressure 
field. 

If the 11 glance” option of calculating the nacelle 
pressure field is used, this subroutine calculates 
what part of a pressure signature from a selected 
nacelle acting at a given sparwise station should 
be deleted. The approximation is made that the 
pressure signature propagates along Mach lines for 
this calculation. 

If deletion of part of the signature is required, 
subroutine NCUT0FF is called. 

CALL NGLANC (Y, YNAC, BETA) 

Y Y station on wing 

YNAC Y station of nacelle generating 

pressure signature 

BETA Mach number parameter, V m 2 - 1 


NCUT0FF 


Subroutine NMERG 

To create a single composite pressure signature 
from multiple superimposed signatures. 

The nacelle pressure field in NACL at a given Y 
station is built up from contributions from all 
nacelles. The pressure signature from a single 
nacelle is contained in arrays BT0 and CT0 (x and 
pressure coefficient) , NW values of each. The 
composite signature becomes XVAL and CPVAL, with 
JNEXT values of each. 

All values of the new pressure signati re 
retained, and are merged with any previous values 
in XVAL and CPVAL, except that the merged pressure 
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USE: 

INPUT: 

SUBROUTINES 

CALLED: 


PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


signature is cut off aft of an input X value (XTE) . 
In addition, the location of the bow shock of all 
pressure signatures in CPVAL is saved (in TSX, 
NWHAT values ) - 

CALL NMERG (JNEXT, NWHAT, XTE) 

JNEXT Number of values in composite 

CPVAL array 

NWHAT Number of combined pressure 

signatures 

XTE Aft-most X value of interest 


TBLU1 


Subroutine NACSET 

To separate nacelles into those above or below 
wing, write nacelle geometry, and calculate nacelle 
wave drag, wetted area, and F(y) function, 

NACSET is used for general bockkeeping involving 
the nacelles, and is called by NACPF. It scans the 
nacelle origins and separates them into those below 
and above the wing. It prints the nacelle input, 
and then zeroes the tables used to sum the nacelle 
interference drag terms. It then calls subroutine 
NB0DY to calculate the nacelle F(y) function, wave 
drag, and wetted areas. 

If all nacelles have the same goemetry, NACSET 
calculates the nacelle F(y) only once, and shifts 
it appropriately for other nacelle origins. 

CALL NACSET 

Nacelle data in common blocks 


NB0DY 


Subroutine N0RDR 

To arrange an arbitrary set of numbers into a 
monotonically increasing array. 

The nacelle bow shock locations contained in TSX 
(NWHAT values) are in random order. For 
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use: 

INPUT; 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


interpolation by TBLUI , they must be in 

monotonically increasing order, and cannot be 
double- valued. N0RDR is used to rework TSX as 
required. 

CALL N0RDR (TSX* NWHAT) 

TSX Array of bow shock X locations 

NWHAT Number of values in TSX array 


None 


Subroutine PINGRT 

To integrate pressure field over surface of 
nacelle. 

Given a pressure signa-ture defined by the 
XVAL, CPVAL arrays, PINGRT is used to inte- 
grate the signature as a buoyancy field over 
a nacelle area distribution. 

CALL PINGRT (J, F0RCE, JNEXT) 

J Index value of nacelle origin 

F0RCE Resultant pressure drag term, 

fc dA 
p 

JNEXT Number of values in XVAL or CPVAL 

array 


TBLUI 


Subroutine NYSET 


To set up semi-span Y 
nacelle pressure field. 

values 

for 

calculation 

of 

A set of Y stations 

on 

the 

right 

hand wing 

are 


identified in NYSET for use in defining the nacelle 
pressure field. This set of Y stations consists of 
the semi- span Y values used in the wing thickness 
pressure calculations (contained in TYB2) plus 
extra stations immediately inbcard and outboard of 
each nacelle centerline. 

The set of Y values is arranged in monotonically 
increasing order by calling subroutine N0RDR. The 
resulting array is TYP, with NYF values. 
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USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


CALL NYSET 

Configuration genunetry passed in through common 
blocks 

N0RDR 


Program L0AE6 

To read input geometry for near- field wave drag 
module. 

L0AD6 reads the input set up by the near- field 
program interface in the geometry module. 

CALL 0VERLAY (NFWD, 6, 2 , 0) 

See input instructions in user’s manual. 


TBLUI 


Program EMPNN 


To calculate 

coefficients for 

thickness 
canard (s) 

pressures 
or fin(s) . 

and 

drag 

EMPNN selects 

empennage 

surfaces in 

turn 

and 


computes thickness pressures and drag coefficients 
utilizing the same logic as wing solution (program 
P916), except interference terms with remainder of 
configuration are neglected. 

CALL OVERLAY (NFWD, 6, 5 r 0) 

Empennage geometry data and basic configuration 
data contained in common blocks. 


GRID, SETGRID, DZEMP, NFWDR, TBLUI, PWRIT 


Subroutine DZEMP 

To calculate surface slopes of canard or fin. 

DZEMP calculates surface slopes of canard or fin 
for all grid elements defining the surface. Slopes 
of grid elements inboard of the side-of-body 
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USES 


INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


station are set to zero. Method parallels that of 
subroutine DZCALC for wing. 

Call DZEMP 

Configuration geometry in common blocks 
GRID, TBLU1 

Subroutine GRID 

To define empennage grid system data. 

GRID is a utility routine, to define grid element 
data at a selected spanwise element row. 

CALL GRID (JBY, XLE, XTE, LLE, LTE, ALE, ATE) 

See description under subroutine WGRID 

None 


Subroutine SETGRID 

To set up grid system for selected canard or fin. 

Given an identified empennage part (defined in 
common block /EPLAN/) , SETGRID sets up the 
corresponding grid element system. Method 

parallels that of routine XFM for wing. 

CALL SETGRID 

Configuration data in common blocks. 


TBLU1 , GRID 


Subroutine PWRIT 

To write out pressure distribution of canard or 
fin. 

Common block /DINT/ contains thickness pressure 
array (ECPE) of canard or fin. PWRIT is called to 
output the data. 
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USE; 

INPUT: 

SUBROUTINE 

CALLED; 


CALL PWRIT 

Configuration and pressure data in common blocks. 
None 


PURPOSE: 

METHOD: 

USE; 

INPUT: 

SUBROUTINES 

CALLED: 


Program WRIT 
To write out input data 

WRIT is used to write the input and pertinent program 
data onto the output file* 

CALL OVERLAY (NFWD, 6, 6, 0) 

Configuration data in common blccks. 

None. 
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2.5 WING DESIGN AND OPTIMIZATION MODULE 

The wing design and optimization module is priirary overlay 8. It 
contains secondary overlays to calculate: 

• The flat wing solution for the input planform 
and Mach number 

• Carry-over pressure distribution on wing 
inboard of side-of-fuselage station- 

• The wing shape and force characteristics for a 
given loading 

• The optimum combination of loadings for least 
drag, subject to imposed constraints. 

A schematic of the wing design and optimization module is 
shown in figure 2.5-1. Major logic of the optimization coding 
is shown in figure 2.5-2. 


PURPOSE: 


METHOD: 

USE: 

INPUT: 


Program TEA - - 2 53 

To call the calculation overlays of the 
wing camber design program and to call the 
interactive graphics displays (if used) . 

Standard FORTRAN statements 

Call OVERLAY (A3 89, 8, 0) 

Input is in DATA OVERLAY (WDEZ , 10, 1) 


Subroutine F0LX0W 

PURPOSE: To provide diagnostic information which identifies 

the entry to and departure from each secondary 
overlay of the wing camber design program. 

METHOD: Standard FORTRAN statements. 


USE: Call FOLLOW (I, J) 

INPUT: The parameters I and J in the calling statement. 

I is the secondary overlay number. If J is 1, then 
overlay I is being entered; if J is 2, then overlay 
I is being departed. 
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FIGURE 2.5-1 -WING DESIGN AND OPTIMIZATION MODULE 
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FIGURE 2.5-2-MAJOR LOGIC OF OPTIMIZATION CODING IN WING DESIGN MODULE 


51 


































Function CP ARB 


PURPOSE: 


METHOD: 


USE: 

INPUT: 


PRINCIPAL 

VARIABLES: 


To define lifting pressure coefficient in the 
arbitrarily defined plant orm region for component 
loading 11. 

At a given span station, linear interpolation is 
used to establish the chordwise coordinate of the 
arbitrary region's leading edge. Then, Cp is 
proportional to distance aft of this location. 

CP = CPARB (XF, YF, CHORD) 

XF the chordwise and spanwise location 

YF on the wing planform in fractions 

of local chord and seirispan at 
which pressure coefficient is to be 
provided. 

CHORD local wing chord at span station YF. 

In addition, a definition of the arbitrary region 
is passed through COMMON/BLOCK 4/. 


CPARB the lifting pressure coefficient 

in the arbitrarily defined region. 
XARB chordwise and spanwise coordinates 

YARB of points defining the arbitrary 

region. 


PURPOSE: 


METHOD: 


USE: 

INPUT: 


Subroutine CPDEF 

To define lifting pressure coefficient as a 
function of chordwise and spanwise planform 
location for each of the seventeen component 
loadings; further, to also define this pressure 
coefficient for an optimum combination of loadings. 

Lifting pressure coefficients fcr loadings 1-11 are 
defined by analytical expressions and for loadings 
12-17 by linear interpolation in tables as a 
function of chordwise and spanwise planform 
location. If a wing-body configuration is being 
designed, then lifting pressure coefficient with- 
in the body is obtained by interpolation within the 
body carry-over tables as a function of planform 
location (using function CARRY0). 

CALL CPDEF (DX, BET AY, CHORD, KCPT, IFLAG) 

DX is distance aft of leading edge in 
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PRINCIPAL 

VARIAELES: 


program units. 

BETAY is the pr oduct of the Mach number 

factor ( J M 2 - 1) and spanwise distance 
from the centerline in program units. 

CHORD is wing chord in program units at 

span station BETAY. 

KOPT is the index for the loading number 

table LOADNO. 

I FLAG is a parameter which offers the 

option of deleting the contribution 
to lifting Cp of the three configu- 
ration-dependent loadings. This 
option is used during camber surface 
calculations which use lifting pres- 
sures, but not during calculation of 
wing surface pressures. 

Input is also provided through the use of 

common blocks. 


AI 

CPBODL 

CPBUPW 

CPDEF 

FLOAD 

LOADNO 

NLOADS 

PCAFRY 

XBODUP 

YBODUP 


a table of loading factors computed 
by program OPTIMUM in OVERLAY 
(WDEZ, 10, 5), one for each of the 
component loadings, that defines the 
optimum combination of loadings, 
a table of lifting pressure coef- 
ficient acting on the wing and 
created by nonsymmetric body 
volume distribution, 
a table of lifting pressure coef- 
ficients acting on the wing and 
created by the body upwash flow 
field. 

lifting pressure coefficient 
(lower surface Cp - upper surface 
Cp) . 

a table of scaling factors set in 
program PI 51 which produces com- 
ponent loading lift coefficients 
of order 1.0. 

a table of component loading 
numbers with index KOPT that are 
input data in program DATA, 
the number of component loadings 
a table of lifting pressure coef- 
ficient in the body carry-over 
region computed by C01IFT. 
the chordwise and s panwise tables 
corresponding to the body upwash 
loading table CPBUPW, in units of 
percent of lccal chord and percent 
of local span, respectively. 
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SUBROUTINES 

CALLED: 


PURPOSE: 

METHOD: 

USE: 

I NFUT : 


PRINCIPAL 

VARIABLES: 


PURPOSE: 


XCPBOD the chordwise and spanwise tables 
YCPBOD corresponding to the body buoyancy 
loading CPBODL r in units of percent 
of local chord and percent of 
local span, respectively. 

YS0B span station of the side of body. 


CPARB, CPINTR, CPNACI, CARRY0 


Function CPINTR 

This program interpolates linearly in chord- 
wise and spanwise directions fcr pressure 
distributions not defined in COMMON blocks. 

Chordwise linear interpolation is used first 
for two span stations bracketing the span 
station of interest; then, linear spanwise 
interpolation is used. 

CP = CPINTR (XF, YF, X, Y, CP, NX, NY, 

NXMAX, NYMAX) 

XF the chordwise and spanwise coordi- 

YF nates of the point on the plan- 

form for which pressure coefficient 
is to be found, in fractions of 
local chord and of serrispan, re- 
spectively. 

X tables of chordwise and spanwise 

Y locations at which pressure coef- 

ficient CP is defined. NX and NY 
values, respectively, are defined 
and NXMAX and NYMAX, respectively, 
are the maximum values of NX and NY. 
CP a two-dimensional array of pressure 

coefficients . 


CPINTR interpolated value of pressure 
coefficient. 


Function CPNACI 

To interpolate linearly both chordwise and 
spanwise for lift pressures due to nacelles 
that are defined in COMMON. The chordwise 
locations at which nacelle pressures are 
defined vary from span station to span 
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station 


METHOD: The span stations bracketing the desired 

span station are located first. Linear in- 
terpolation chordwise is completed at these 
two span stations for the desired chordwise 
location, and finally, spanwise interpolation 
is carried out at the desired chordwise 
location. 

USE: CP = CPNACI (XF, YF) 

INPUT: XF chordwise and spanwise locations 

YF on the wing planform in fractions 

of local chord and seirispan, re- 
spectively. 

In addition, the nacelle pressure field CP 
and the spanwise and chordwise locations at 
which it is defined, Y and X, are defined by 
COMMON NPF. 

PRINCIPAL 

VARIABLES: CP a two-dimensional table of nacelle 

pressure field defined in the 
analysis program. 

CPNACI interpolated value of the nacelle 
pressure coefficient. 

X a two-dimensional table of chord- 

wise locations corresponding on a 
one-to-one basis with the CP table. 
Y a one-dimensional table of semispan 

stations corresponding to the span- 
wise variations of both X and CP. 


Function CPTDEF 

PURPOSE: To interpolate linearly both chordwise and 

spanwise for pressure coefficient due to wing 
thickness and defined in COMMON. 

METHOD: Linear interpolation is carried out chordwise 

at the two span stations bracketing the de- 
sired span station, and is followed by 
linear spanwise interpolation. 

USE: CPT = CPTDEF (XF, YF) 

INPUT: XF chordwise and spanwise stations on 

YF the wing planform in fractions of 

local chord and semispan, respect- 
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ively. 


PRINCIPAL 

VARIABLES: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

OUTPUT: 

SUBROUTINE 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 


In addition, the wing thickness pressures 
are defined in a two-dimensional array CPT 
in COMMON block PLIM; the chordwise locations 
X and spanwise locations Y corresponding to 
the elements in CPT are also defined in 
COMMON. 


All defined above. 


Subroutine JXTIME 

To print out time since job began execution 
and to print out the time increment since 
the immediately preceding call to this 
subroutine. 

Interrogate the system timing subroutine 
SECOND. 

CALL JETIME 

None 

Time increment and time mentioned above, 
both in decimal seconds. 


SECOND 


Function CARKY0 

To compute lifting pressure coefficient within the 
body carry-over region by interpolation within the 
tabular data provided by C0LIFT as a function of 
planform location. 

Linear interpolation chordwise and spanwise. 

CP = CARRY0 (DX, BY, CH0RD, K0PT) 

DX is distance aft of the leading edge in 

program units. 

BY is the product of the Prandtl- 

Glauert Mach number parameter and 
the spanwise distance from the 
centerline, in program units. 

CH0RD. is wing chord in program units at span 
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K0PT 


PURPOSE: 

METHOD: 

USE: 

INPUT: 


PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 


PURPOSE: 
METHOD : 


station BY. 

is the index for the loading number in 
table L0ADN0. 

Input is also provided by the use of common blocks. 


Subroutine WDGRID 
To define wing grid system. 

WDGRID is a utility routine to identify the grid 
elements defining, the wing planform at a selected 
spanwise element row. 

CALL WDGRID (JBY, XLE, XTE , LLE, LTE, ALE, ATE) 

JBY selected spanwise element row 

XLE wing leading edge coordinate at JBY 

XTE wing trailing edge coordinate 

LLE wing leading edge element 

LTE wing trailing edge element 

ALE element fraction at leading edge 

ATE element fraction at trailing edge 


Program DATA 

To read input data for the wing design module. 

Program DATA reads the input data (set up by the 
wing design interface in the geometry module) , and 
writes it out according to the print code set in 
the input. The wing grid system is also 
established (in subroutine TRNSFM) . 

CALL 0VERLAY (WDEZ, 8, 1) 

See User's Manual 


TRNSFM, MAXMIN, REGRXD, JLTIME 


Subroutine MAXMIN 

To identify and print both the maximum and 
the minimum elements in a two-dimensional 
array. 

Standard FORTRAN library subroutines. 
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USE: 
INPUT 2 


OUTPUT: 


CALL MAXMIN (A, NX, NY, NX MAX, NYMAX) 

A a two-dimensional array with maxi- 

mum dimensions NXMAX and NYMAX 
containing (NX) (NY) values to be 
searched for maximum and minimum, 

THEMAX The maximum element of A 

THEMIN The minimum element of A. 


PURPOSE: 

METHOD: 


USE: 


INPUT: 


PURPOSE: 

METHOD: 


Subroutine REGRID 

To define grid system data when the RESTART 
option is used. 

When the RESTART option is used, all basic 
loading data is input from cards (or tape) 
and grid system calculations normally made 
(which will be needed) are bypassed. 

REGRID sets up the appropriate arrays, which 
include the wing chord at each of the span- 
wise calculation stations, and the per cent 
chord and per cent wing length associated 
with each grid element along these spanwise 
stations. 

CALL REGRID 

COMMON block BLOCK 1 is used to input a 
definition of the wing planform and the 
associated parameters required to define 
the Mach box grid system. 

Subroutine TRNSFM 

To convert input data to program units and set up 
wing grid system. 

The wing is represented in the program by a set of 
rectilinear elements, with the number of semi-span 
element rows given by input TN0N. TRNSFM 
interpolates the planform at the centerline of each 
element row to define the leading edge and trailing 
edge values, converts them to program scale (using 
RATIO = Ywing tip*BETA/TN0N) , and stores them in 
arrays TXLE and TXTE. Special values XLE0 and XTE0 
define the wing centerline leading edge and 
trailing edge. 
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USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 


If the parabolic apex option is selected in the 
input (YSN00T>O.) , the wing leading edge out to 
YSN00T is altered to a parabolic shape. 

If a fuselage is input with side-of-f uselage 
station greater than 0., the grid system element 
row (NS0B) corresponding to the side-of-f uselage 
station is identified. TRNSFM then checks to see 
that NS0B is included in the spanwise camberline 
row array (TJBYS) , adding it if necessary. 

CALL TRNSFM 

A definition of the wing planform in physical units 
is passed to TRNSFM by COMMON block BLOCK 1 and 
block SNOOT. 


TELU1 


Program PI 51 

To calculate the flat wing leading; and to set 
normalizing factors for component loadings 1-10 so 
that their lift coefficients will be approximately 

1 . 0 . 

The method of the flat wing solution is essentially 
the same as described for the lift analysis module 
(using the equations given in the theory document, 
part 1) , except that only the wing is present (no 
fuselage, nacelles, etc) . The wing is scanned from 
front to back and centerline to right hand wing 
tip, computing the pressure coefficients (CP) for 
all field point elements (LSTAR, NSTAR) . The 9- 
point smoothing equation is applied after all 
pressure coefficients are calculated. 

The lifting pressure distribution is calculated 
over the surface of the wing to obtain lift 
coefficient (SCL9) , drag coefficient (SCD9) , 
dCm/dCL (DCMCL) and drag-due- to-li ft factor (KF) , 
based on input reference geometry. 

For the given planform, normalizing factors (array 
FL0AD) are then calculated which will produce a 
lift coefficient of approximately unity when used 
with each of the analytically defined basic 
loadings (1-10) in program P91615. 

CALL 0VERLAY (WDEZ, 8, 2, 0) 
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INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 


SUBROUTINES 

CALLED: 


Configuration data in common blocks 


JLTIME 


Program C0LIFT 

To calculate carry-over lift distribution in fuse- 
lage region. 

C0LIFT creates a definition of the carry-over lift 
distribution for all selected loadings, to be 
applied to the wing planform region covered by the 
fuselage (if there is one) . The calculation 
procedure utilizes the analysis form of the lifting 
pressure solution (described in the lifting 
analysis program writeup) for spanwise element rows 
inside the side-of-body station. The resultant 
carry-over lift distributions are stored in array 
PCARRY. Carry-over lift distributions are not 
computed for the configuration dependent loadings 
(loadings 15, 16, and 17) since those are already 
of the carry-over type. 

C0LIFT also integrates lift distributions to get 
drag, lift, and pitching moment of loadings acting 
on fuselage, which are stored in common block 
/FUSCX/ . 


CALL 0VERLAY (WDEZ, 8, 3,0) 


NS0B 

TXB, TZB 
B0DCL (I) 

B0DCM (I) 
B0DCM (I) 

PCARRY 

(I,M,NYS) 


USEB0D 


side-of-body element row 

fuselage camber shape 

carry-over lift of I-th loading act 

ing on fuselage. 

carry-over drag of I-th loading 

carry-over pitching moment of Ith 

loading 

carry-over pressure distribution 
due to I-th loading. NYS spanwise 
rows given in array KYS, M chord- 
wise percentages given in array 
SXPC. 

logical flag (carry-over lift 
calculated if USEB0D = TRUE) 


CPDEF, JLTIME, TBLU1, VSDGRID 
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OVERLAY (WDEZ, 10 , 4) 
Program P91615 

PURPOSE: 


METHOD: 


QE^ion Q)_: 

In the calculation of the camber surface and force coefficients, 
a series of semi-span stations are selected in the program input 
(TJBYS) . The program then picks a loading, calculates the 
required surface slopes for all elements at each TJBYS station, 
and integrates the slope distribution to obtain the camberline. 
In these calculations, each element along TJEYS is identified by 
the nomenclature LSTAR (x) and NSTAR (3y) , as described in the 
theory document. Z ordinates at specified points on the planform 
are interpolated and saved if Z constraints are to be applied. 

After the camberlines and sectional force coefficients of all 
spanwise stations have been computed, the characteristics are 
integrated spanwise to obtain lift coefficient (CL) , drag 
coefficient (CD) , and pitching moment coefficient (CMAPEX) . These 
are converted to input reference geometry basis and become CLR and 
CDR, lift and drag coefficients. The pitching moment coefficient 
is adjusted to the value at zero lift, and becomes CM0R. Drag- 
due-to-lift factor is labeled KE = CDR/ (CLR) 2 . The interference 
drag coefficients are stored in array CDI. 

After all calculations are completed for a given loading, the 
force coefficients are stored in array TDRAG, and the process 
repeated until all loadings have been used. 


To calculate the aerodynamic characteristics of a 
specified lift loading and the camber surface 
required to support it. Both component loadings 
and combinations of component leadings are handled. 
If requested, all data for the RESTART option are 
punched in this program. 

Program P91615 solves for the camber surface 
required to support a specified loading and the 
associated force coefficients, using the equations 
given in the theory document (part 1) . The program 
is actually used in two ways: 

1) To calculate the force coefficients and 
interference drag characteristics of a set of 
basic loadings. 

2) To calculate the camber surface for an optimum 
combination of loadings. 

The program code used to define the usage of P91615 
is K0PTI. If K0PTI is greater than the number of 
basic loadings used, then option (2) above is 
employed. 
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The force coefficients and interference drag coefficients of all 
loadings are then converted to the component and interference 
forms used in the matrix solution described in the theory 
document, and stored in common blocks to be passed to program 
0PTIMUM. 

Finally, the RESTART data are written onto tape 3, and also 
punched into cards, if requested. These data consist of all 
component and interference drag terms, the configuration- dependent 
loadings {if used), and grid system data calculated by P91615. 


QE£igr?_i21: 

If option 2 was selected, involving the calculation of the wing 
shape for an optimum combination of loadings, the calculation 
sequence is the same as if a basic loading was being used. 
However, the interference characteristics are not required. The 
resulting camber surface is stored in common block/CAMBER/ in the 
following form: 


TFCT 

PYB2 

TZ0RD (x,y) 
IFZC 


percent chord array, N0PCT values, 
semi-span array, JBYMAX values, 
camber surface 2/C values, in per cent. 
1 , to denote 2/C in per cent. 


The camber surface will also be punched into cards, if requested 
in the program input. 

USE: CALL OVERLAY (WDEZ, 8, 4) 

INPUT: All input is handled by COMMON blocks, which 

pass the required data from DATA. 

SUBROUTINES 

CALLED: CDWONN, CPDEF, CPINTR, CPNACI , CPTDEF, 

HEADER, JLTIME, SM00TH3, TBLU1 


PURPOSE: 

METHOD: 


USE: 


Function CDWONN 

Tc calculate the axial force acting on 
nacelles due to the wing-lift flow-field. 

The wing lower surface lifting pressures are 
projected downward in a vertical plane and 
aft along Mach lines to the nacelle locations 
where their product with nacelle incremental 
frontal area is numerically integrated to 
produce a nacelle axial force. 

CDN = CDWONN (K0PT) 
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INPUT: 


K0PT 


component loading index 


PRINCIPAL 

VARIABLES: 


SUBROUTINE 

CALLED: 


PURPOSE: 

METHOD: 

USE: 

INPUT: 


In addition, wing planform information is 
passed to CDWONN by COMMON block BLOCK1 
and nacelle geometry is passed by block 
block 1 1 . 


CPTERM 

XLEN 

XTEN 

CRDN 

RNAC 

XNAC 


XYZ 


wing lower surface pressure coef- 
ficient due to wing lift only, 
longitudinal coordinates of the 
wing leading and trailing edges at 
the nacelle span station, and the 
corresponding wing chord, 
two-dimensional arrays defining 
nacelle radius as a function of 
nacelle longitudinal station. 

The first parameter is the nacelle 
number index. 

a two-dimensional array specifying 
the coordinates of each of the 
nacelles* forward-most point. The 
first parameter is the nacelle 
number index, and the second para- 
meter defines the coordinate being 
referenced - i.e. , first value is X, 
second is Y, and third iz Z. 


CPDEF 

Subroutine HEADER 

To write out loading number and case 
identification. 

Standard FORTRAN statements. 

CALL HEADER (K0PT) 

K0PT component loading index. 

In addition, the loading number index, 
loading name, and case identification are 
passed via COMMON statements. 

Titling information for each component 
loading prior to its aerodynamic and camber 
analysis. 


OUTPUT: 



Subroutine SMOOTH 3 


PURPOSE: 

METHOD: 


USE: 

INPUT: 


PURPOSE: 

METHOD: 


To apply 3-point smoothing to a selected 
range of elements within a one-dimensional 
array . 


Each element to be smoothed is replaced 
according to the following algorithum: 

i (y i-i 2y i + y i+i' 

where y i-l and y i+l are the immediate 
neighbors of y i before smoothing. If 
is the first element in the array to be 
smoothed, then 

y i = i <2y i + y i+ i> 


and if y- x is the last. 


y i =i (y i-i + 2y i> 


CALL SMOOTH3 (A, IFIRST, ILAST, N, NARRAY) 


A 

IFIRST 

ILAST 

N 

NARRAY 


the one-dimensional array of data 
to be smoothed. 

first and last elements in A to be 
smoothed. All elements between 
will be smoothed. 

the number of times the smoothing 

algorithum is to be applied. 

the maximum number of elements in A. 


OVERLAY (WDEZ, 10, 5) 

Program OPTIMUM 

To define various optimum combinations of lift 
loadings in terms of their lead strength factors 


Lagrange* s method of undermined multipliers {as 
described in Part Is Theory) , as a function of the 
aerodynamic characteristics of each of the 
component loadings and their mutual interferences. 
In program OPTIMUM, the DO loop on statement 730 
(index ILOOP) is used twice only if a constraint on 
pitching moment coefficient is used. For ILOOP = 
1, program OPTIMUM produces a solution for minimum 
drag with only a lift constraint (this solution 
also includes constraints on ordinates if they are 
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used); 21 solutions of the drag-due-to- lif t factor 
(KE) and zero lift pitching moment coefficient 
(CMO) defining the design "bucket” plot; and if 
requested, a solution for minimum drag with lift 
coefficient and wing upper surface pressure 
constraints (this solution also includes 
constraints on ordinates if they are used) . If the 
wing pressure criteria are satisfied by the first 
solution, the latter solution is set equal to the 
first. For ILOOP = 2, program OPTIMUM adds a 
constraint on pitching moment coefficient at zero 
lift to both the solution with lift coefficient 
constraint and the solution with lift and pressure 
constraints . 

Within the ILOOP loop, the left-hand-side solution 
matrix AMAT and the right-hand- side solution matrix 
BMAT are calculated first, corresponding to figure 
4.4-4 of the theory document (part I). The left- 
hand-side matrix is stored in ATEMP for 
multiplication with the solution as a test of its 
accuracy; this multiplication should produce the 
right-hand-side. 

Subroutine MATINV is called tc solve for the load 
strength factors A^ and the Lagrange multipliers, 
corresponding to the left-hand-side column matrix 
in figure 3.4-3 of the theory document- The 
solution load strength factors Ap are stored in two 
arrays - the AI array so that the current solution 
is defined in the wing lifting pressure subroutine 
CPDEF, and in the TAI array, so that the solution 
will be defined for program P91615 after program 
OPTIMUM has been exited. The array TAI has 
capacity for four sets of load strength factors 
(under IDUM) , corresponding to the four types of 
available solutions as follows: 


IDUM Lift 

Constraint 


Pitching 

Moment 

Constraint 


Wing 
Pressure 
Constraint (s) 


1 

yes 

no 

no 

2 

yes 

yes 

no 

3 

yes 

no 

yes 

4 

yes 

yes 

yes 
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Lift coefficient CLSOL, drag coefficient CDSOL, and 
pitching moment coefficient CMOS are computed in 
program units from the load strength factors and 
the aerodynamic coefficients of the component 
loadings. These parameters are then converted to 
the input reference geometry basis, CLR, CDR, and 
CMOR, respectively. Values of KE and CMO from 
previous solutions (if any) for IDUM = 4 are 

shifted one location toward the rear of the 
"bucket” plot arrays CKE and CM2ERO, and the 
current solution data are stored in these arrays. 

The wing upper surface pressure coefficient CPUS 
and the corresponding longitudinal gradient DCPDX 
are calculated next and compared everywhere with 
the user defined limiting values, CPLIM and CPGLM, 
respectively. The minimum difference of CPUS 
CPLIM is identified (CPMIN) and the maximum 

difference of DCPDX - CPGLM is identified (GRADM) , 
along with their respective planform locations. 
This completes the solution corresponding to IDUM 
- 1 above. 

Then, if ILOOP = 1 and if no pressure constraints 
have been applied, values cf KE and CMO are 
generated by the DO loop on 630 for the "bucket" 
plot. This solution parallels the one just 

described, except that a constraint on design 
pitching moment coefficient at zero lift CMOD is 
added and is varied through 21 values. The range 
of values of CMOD depends on design lift 

coefficient CLD2IN, and is centered on the pitching 
moment coefficient corresponding to the solution 
for IDUM = 1 , if available, or on zero- Values of 
CMOD are truncated for plotting convenience. 

After the bucket plot data are generated and 
stored* program OPTIMUM tests tc see if pressure 
constraints on the wing upper surface have been 
requested and whether they are necessary, if 
requested. If both tests are positive, then a loop 
to statement 10 is used to apply constraints on 
wing pressure. This loop is within the loop on 730 
(index = ILOOP) • The loop which tests wing upper 
surface pressure functions as follows. First, the 
critical pressure gradient is tested to see if it 
is satisfactory. If it is rot, a constraint on 
pressure gradient is applied at the planform 
location where the gradient criterion is violated 
by the greatest margin. At this location, the 
gradient is constrained to 75 percent of the 
gradient criterion; this provides a margin to 



prevent the gradient from exceeding the limit in 
the immediate vicinity of the constraint. The 
program then loops to statement 10 and the process 
is repeated until the overall gradient criterion is 
satisfied. 

The program then tests critical pressure level 
against its acceptable level. If it is not 
acceptable, a pressure level constraint is added at 
the planform location where the pressure level 
criterion is most strongly violated. The lifting 
pressure coefficient at the critical planform 
location is calculated which just meets the upper 
surface level criterion, and the lifting pressure 
is constrained to 95 percent of this value. The 95 
percent factor provides a margin to prevent the 
pressure from exceeding the criterion in the 
immediate vicinity of the constraint location. The 
program then loops to statement 10, and a new 
solution is generated. 

It has been found that both pressure gradient and 
pressure level constraints imposed early in the 
design cycle can become unnecessary, and even 
undesirable, when later constraints are imposed. 
This situation can be detected in the solution by 
noting that the algebraic signs of the pressure- 
constraint Lagrange multipliers change from 
positive to negative when the constraints become 
unnecessary. 

The design program takes advantage of this by 
testing the Lagrange multipliers of the pressure 
constraints for each solution cycle. If there are 
unnecessary constraints, they are removed, and the 
solution cycles back to statement 10 with fewer 
pressure constraints. 

In summary, then, each solution is tested for wing 
upper surface pressure gradient and pressure level, 
and for negative pressure constraint Lagrange 
multipliers. The hierarchy is; 

(1) Test for unnecessary pressure constraints 

(2) Test most critical pressure gradient. 

(3) Test most critical pressure level. 

If one of these is violated, appropriate action is 
taken, and the program loops back to statement 10. 
Throughout this cyclic process, each solution 
includes constraints on 2 ordinates if such 
constraints have been requested. 



It is certainly possible to define wing upper 
surface pressure criteria which cannot be 
satisfied, either through input errors or simply 
through ignorance of, say, wing thickness pressures 
on the first run of a new configuration in the 
system. Should this occur, the program will add 
pressure gradient and/or pressure level constraints 
until it reaches one of two possible constraint 
maxima. 

One maximum on the number of constraints is imposed 
by the number of loadings used. The total number 
of constraints, including those on lift, pitching 
moment, configuration-dependent loadings, 
ordinates, and pressure, can at most be equal to 
the number of loadings. This situation is 
undesirable for it leaves no degrees of freedom for 
drag minimization; consequently, this limit has 
been set so that two degrees of freedom remain free 
for drag minimization (for small numbers of 
loading, this is reduced to one). 

The second maximum imposed on the number of 
pressure constraints is dictated by the number of 
loadings that are free to influence longitudinal 
pressure gradient. A maximum of ten loadings do so 
- loadings 2, 5-9, and 11-14. It has been thought 
desirable to leave one degree of freedom for drag 
minimization for gradient constraints. For this 
purpose, the number of permissible gradient 
constraints is further reduced by one whenever a 
constraint is imposed on C rao since C mo constraints 
are satisfied primarily by the same x-dependent 
loadings used to satisfy gradient criteria. 

If the program reaches a pressure constraint 
maximum, it checks to see if gradient constraints 
have been imposed. If one or more gradient 
constraints have been imposed, the program 
arbitrarily increases the gradient criterion table 
by 20 percent, and begins anew with no pressure 
constraints. This process is also cyclic and can 
be repeated up to 50 times before halting with the 
solution produced by the last cycle. 

If no pressure gradient constraints have been used, 
the program halts immediately upon reaching one of 
the two loading maxima. 


For I LOOP = 2, a pitching moment constraint is 

added to the solution, and lift coefficient. 



USE: 

INPUT: 

PRINCIPAL 

VARIABLES: 


ordinate and pressure limiting constraints function 
as they did for ILOOP = 1. The bucket plot 

calculations are omitted. 

Major logic of the optimization sequence in the 
wing design module is illustrated in Figure 2.5-2, 

CALL OVERLAY (WDEZ, 8, 5) 

All input is by means of COMMON blocks. 


AI 

AMAT 

EMAT 

CDIJ 

CDSOL 

CDWON 

CLDSIN 

CLI 

CPBODL 


CPBODU 

CPBUPW 

CPLIMIT 

CPNAC 

GRADL 


GTEMP 


MAXN0G 


NCMAX 

NCPCON 


load strength factor Aj. 
left-hand-side solution matrix 
right-hand-side solution matrix 
interference drag coefficient 
solution drag coefficient 
interference drag coefficient of 
wing lift acting on nacelles, 
design lift coefficient 
the ith component lift coefficient 
lifting pressure coefficient due 
to unsymmetric body volume 
distribution 

body pressure coefficient acting 
on the wing upper surface 
lifting pressure coefficient due 
to the body upwash loading 
wing upper surface limit pressure 
coefficient (on input table) 
wing lifting pressures due to 
nacelles 

the array of maximum longitudinal 
pressure gradient on the wing upper 
surface actually used in the wing 
design module. This array is 
increased by 20 percent for 
those solutions having maximum 
number of pressure constraints, 
including at least one on 
gradient. 

the array of user-specified maxi- 
mum pressure gradient. This array 
is copied into GRADL on entry into 
OPTIMUM. 

the number of loadings that contri- 
bute to longitudinal pressure grad- 
ient and that are free to vary, 
maximum number of solution con- 
straints 

number of pressure coefficient 
constraints 
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OUTPUT: 


NGC0N 

NLCON 

NXCMAX 

N2CON 

OPTION 

TAI 

TBETAY 

TDXM 

CHORDT 

CPLMX 

USEBOY 

USEBUP 

USECMC 

USECPL 

USEOPT 

USENAC 

US ESC 

X2CON 

Y2CON 

SCON 

21(17,5) 


the number of solution constraints on wing 
upper surface longitudinal gradient- 
number of direct solution loading 
constraints due to use of configu- 
ration-dependent loadings 
the maximum number of pressure gradient 
constraints. 

number of ordinate constraints 
integer array controlling extent 
of four types of solutions 
up to four sets of load strength 
factors Ai , corresponding to the 
four types of solutions, 
stored parameters which define a 
solution lifting pressure con- 
straint, namely, spanwise location 
chordwise location, local chord, 
and allowable lifting pressure 
coefficient, respectively, 
a logical flag indicating use of 
body buoyancy loading if true 
a logical flag indicating use of 
body upwash loading if true 
a logical flag indicating use of 
pitching moment constraint if true 
a logical flag indicating use of 
wing upper surface pressure 
limiting if true 

a logical flag indicating use of 
wing thickness pressures if true 
a logical flag indicating use of 
nacelle buoyancy loading if true 
a logical flag indicating use of ordinate 
constraints if true, 
longitudinal and lateral planform 
locations of the ordinate constraints in 
physical units. 

constrained values of ordinate 

corresponding on a one-to-one basis with 
X2CON and Y2CON. 

a two-dimensional array of up to five 
ordinates for each of the seventeen 
component loadings. 


Essential output is the set of up to four solution 
definitions in terms of the loading strength 
factors Ai. These are passed by COMMON blocks back 
to the camber surface calculation overlay, OVERLAY 
(WDEZ, 10, 2). In addition, varying amounts of 
information about the solutions are printed, 
depending on the choice of the print control 
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SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 


SUBROUTINES 

PURPOSE: 

METHOD: 

USE: 

PURPOSE: 


parameter. These are described in more detail for 
the example case in the User*s Manual, part 2. 


CPDEF, CPINTR, CPNACI, CPSL0P, CPTDEF, JLTIME, MATINV, 
SLOPET 


Subroutine CPSL0P 


To provide the longitudinal gradient of wing upper 
surface pressure for component loadings 
individually and when combined in a solution. This 
subroutine parallels subroutine CPDEF in structure. 

Standard FORTRAN statements. 


CPSL0PE = CPSL0P (DX, BETAY, CH0RD, K) 


DX 

BETAY 


CH0RD 

K 


distance aft of the wing leading edge in 
program units. 

the product of the Prandtl-Glauert Mach 
number parameter and the spanwise dis- 
tance from the centerline, in program 
units. 

wing chord in program units at span 
station BETAY. 

index for the loading number table 
L0ADHO. 


Input is also provided by means of COMMON blocks. 


USED: 


CPARB, SLOPEB, SL0PEN, SLOPEU. 


Subroutine MATINV 

To solve a set of linear, simultaneous equations. 
This is a NASA-LRC library subroutine. 

See LRC library. 

CALL MATINV (A # N, B, M, DETERM, IPIVOT, 

INDEX, NMAX, I SC ALE) 

Function SL0PEB 

To provide the longitudinal gradient of the wing 
upper surface pressure coefficient due to the body- 
buoyancy configuration loading. 
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METHOD: 


USE: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 


PURPOSE: 

METHOD: 

USE: 

INPUT: 


PURPOSE: 

METHOD: 

USE: 


Isolation of the body buoyancy contribution to wing 
upper surface pressure coefficient at the chordwise 
stations which bracket the input chord station. 

The gradient is simply the difference of these 
pressures divided by the difference in longitudinal 
distance, all at the proper span station. 

CPB = SL0PEB (XF, YF, C, J) 

Function SL0PEN 

To provide the longitudinal gradient of wing upper 
surface pressure due to the nacelle pressure distri- 
bution when used as a camber loading. 

As for SL0PEB. 

CPN = SL0PEN {XF, YF, C ) 

XF chordwise planform location as a fraction 

of local chord. 

YF spanwise planform location as a fraction 

of semispan- 

C chord at span station YF. 

Tabular input is also provided by the COMMON block NPF. 
Function SL0PET 

To provide the longitudinal gradient of wing thickness 
pressures. 

As for SL0PEB. 

CPT = SL0PET (XF, YF, CH0RD) 

XF chordwise planform location as a fraction of 

local chord. 

YF spanwise planform location as a function of 

semispan. 

CH0RD wing chord at span station YF. 

Wing thickness pressures are provided in tabular form 
by array CPT in COMMON block PIXM. 

Function SE0PEU 

To provide the longitudinal gradient of the wing 
upper surface pressure due to the body upwash con- 
figuration-dependent loading. 

As for SL0PEB. 

CPU ~ SL0PEB (XY, YF, C) 
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I 


INPUT: XF chordwise planform location gives as a 

fraction of local chord. 

YF spanwise planform location as a fraction 

of semispan. 

C wing chord at span station YF- 

The loading due to body upwash is provided in array 

CPBUPW by COMMON block CPBUPW1 - 
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2.6 LIFT ANALYSIS MODULE 

The lift analysis program is primary overlay 7. The program is 
divided into separate elements to read input, transform input into 
program units, and perform the lifting pressure calculations as 
shown schematically in figure 2.6-1. 

The graphics displays are called as shown in figure 2.6-1. 


Program TEA201 

Primary level of lift analysis module 

TEA201 sets up the calculation sequence for the 
drag-due-to-lift analysis program. The calculation 
loops are: 


D0 50 JD0 Mach number loop, repeated for 

each Mach number. 

D0 40 MLIMT Pressure limiting loop. LIMIT 
angles of attack, if limiting 
requested. 

D0 30 JCALP Canard angle of attack loop. 

Repeated for each canard 
alpha, if canard is present. 


CALL 0VERLAY (ANLZ, 7, 0, 0) 
See user’s manual. 


FCALC 


Subroutine DUBINT 

PURPOSE: To perform double interpolation in array 

METHOD: Given a two-dimensional array, DUBINT performs 

double linear interpolation for an answer at a 
specified location in array. 

USE: CALL DUBINT (XI, Y1 , TX, TY, NX, NY, TBL, MX, MY, 

ANS) 

INPUT: XI X location 

Y1 Y location 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 


PURPOSE: 

METHOD: 
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input 


TEA-201 
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FIGURE 2.6-1 .—LIFT ANAL YS/S MODULE SCHEMATIC 






SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 


TX 
TY 
TBL 
MX, MY 
ANS 


X array, NX values 

Y Array, NY values 

Table being interpolated 

Size of TBL in TBL DIMENSION statement 

Interpolated answer at XI, Y1 


None 


Subroutine GETR 

To calculate influence factor for specified grid 
element. 

Given an element located in the forecone from a 
selected field point element, GETR provides the 
corresponding influence factor. (The influence 
factor equation is discussed in the theory 
document, part 1 ) . 

CALL GETR (LSTAR, LVAR, NDIF, E) 

LSTAR Field point element 

LVAR Specified grid element 

NDIF f'NSTAR - N] 

R Influence factor 


None 


Subroutine GGRID 

To define leading edge and trailing edge data for 
grid system for specified semi- span element row. 

Given a wing semi-span station, GGRID defines the 
leading edge and trailing X values, the 
corresponding grid elements, and the associated 
element fractions. GGRID is used for wing, canard 
or horizontal tail. 


CALL GGRID (JBY, XLE, XTE, LLE, LTE, ALE, ATE, 
NK) 


JBY 

XLE, XTE 
LLE, LTE 
ALE, ATE 


Semi-span element row 

Leading and trailing edge X-locations 
of planform at JBY. 

Grid elements corresponding to XLE, 
XTE. 

Fractions of elements defining plan- 
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SUBROUTINES 

CABLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


form at LLE, LTE. 

NK Code to identify configuration component 

1 = Canard 

2 = Wing 

3 = Horizontal tail 


None 

Subroutine RCALC 

To calculate a standard set of influence factors 

The influence factors used in the lift analysis 
program are a function only of the relative 
position of the field point element and the 
influencing element. RCALC is used to precalculate 
a standard set of influence factors for use in the 
program to reduce computer tirre associated with 
repeated calculations. The factors are stored in 
array TRSAVE. 

CALL RCALC 

None 


None 


Program FUSLGE 

To calculate upwash field of fuselage and to 
calculate isolated fuselage force coefficients. 

Fuselage is used to calculate the wing-fuselage 
intersection (subroutine INTSEC) , then to calculate 
the isolated fuselage lift distribution and force 
coefficients using slender body theory, and then to 
calculate the fuselage upwash field acting in the 
plane of the wing, canard, or horizontal tail. 

The equations used in the fuselage lift 
distribution and upwash field calculations are 
given in the theory document (part 1) . The 
fuselage forces in the presence of the wing 
downwash field are later repeated in subroutine 
FUSCF under overlay (7,5). The upwash field of 
canard, wing, or horizontal tail is defined by an 
array of upwash value at specified chord and semi- 
span percentages: 
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CANARD 

WING 

TAIL 

% chord 

TXUPW 

TXPW 

TXUPW 

Semi -span 
fraction 

TYCCW 

TYUPW 

TYHTW 

Upwash angle 

TUPCC 

TUPWC 

TUPHC 

Upwash angle 

TUPCF 

TUPWF 

TUPHF 


per deg. angle 
of attack 


The upwash angle calculations are performed 
in subroutine UPWASH. 


USE: CALL OVERLAY (ANLZ, 7 , 1, 0) 

INPUT: Geometry definition contained in common blocks. 

M1IMT = loop index variable froir (7,0) overlay. 


SUBROUTINES 

CALLED: INTSEC, UPWASH 


PURPOSE: 

METHOD: 


USE: 

INPUT: 


Subroutine INTSEC 

To locate wing-fuselage intersection 

If input SYMM is less than zero, then the analysis 
program is to solve for the wing-fuselage 
intersection in order to define the exposed and 
"carry-over" wing pieces. INTSEC selects each 
percent chord line of the wing camber surface 
definition and locates the intersection between 
wing and fuselage. The fuselage area distribution 
is considered to be made up of circular cross- 
sections in the intersection calculations. The 
resultant intersection is stored in common 
block/WBINT/, with X, Y, and 8 values in arrays WX, 
WY, and W2. 

CALL INTSEC. 

Configuration geometry contained in common 
blocks. 


SUBROUTINES 
CALLED: TBLU1 
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I 


PURPOSE: 

METHOD: 


USE: 

INPUT: 


SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


Subroutine UPWASH 
To calculate fuselage upwash 

UPWASH calculates fuselage upwash angle in the 
plane of the canard, wing, or horizontal tail, 
using the slender body equations discussed in the 
theory document (part 1). Upwash angles are 
computed for a series of percent chord values at 
selected semi-span stations. 

CALL UPWASH (Y, DELX, I, L) 

Configuration geometry contained in common blocks, 
plus: 

Y semi-span y station 

I span storage index for upwash array 

L variable defining component 
1= canard 
2= wing 
3= tail 


TBLU1 


Program WING 

To calculate lifting pressure distribution on wing. 

A schematic of WING is shown in figure 2.6-2. 

The equations used in calculation of the wing 
lifting pressures are given in the theory document 
(part 1). The program scans the wing grid system 
from front to back (D 0 470) and from the centerline 
to right hand wing tip (D0 450) , locating field 
point elements on the wing or canard. When a field 
point element (LSTAR, NSTAR) is located, the 
program (D0 200) computes the upstream influence of 
elements located in the Mach forecone from LSTAR, 
NSTAR. The local pressure coefficient (CP) is then 
computed, with the fuselage upwash added to the 
local surface slope to obtain the effective element 
angle of attack. If the field point element is 
located inside the side-of-fuselage station, the 
element angle of attack is set to zero- (Either of 
the two pressure coefficient smoothing options may 
be used in these calculations, controlled by ->ut 

variable SM0G0.) After CP is calculated, pi e 


79 



Slopes for grid system 


Nacelle pressure signature 


Planform definition (grid system) 


Influence factors 


Limiting pressure check 


Pressure coefficient smoothing 


Lifting pressure tabulation 


Downwash on horizontal tail 


Wing-on-nacelle term 


FIGURE 2.6-2.-PROGRAM WING SCHEMATIC 
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USE: 

INPUT: 

SUBROUTINES 

CALLED: 


limiting is applied if input LIMIT is greater than 
zero, using subroutine CKLIM. 

After all pressure coefficients are calculated, 
they are smoothed (D0 630) and integrated over the 
configuration surface to get lift, drag, and 
pitching moment. At the same time, the local 
nacelle pressure coefficients (CPNAC) and 
asymmetric fuselage volume pressure coefficients 
(CPASYM) are superimposed. Pressure coef f f icients 
are applied to the wing slopes in the exposed wing 
area portion, and to the fuselage slopes interior 
to the side-of-fuselage station. 

Corresponding calculations and summations are 
carried on simultaneously for the flat wing at 1 
degree angle of attack. The interference drag term 
of flat wing pressure coefficients on the cambered 
wing slopes is also computed. Separate summations 
carry the nacelle drag, lift, and pitching moment 
and the configuration streamwise and spanwise lift 
distributions. 

The wing pressure coefficient arrays (TWCP and 
TWCPF) are then interpolated over the wing planform 
for the output pressure summary (subroutine CPTAB) . 

The downwash due to the wing acting in the plane of 
the horizontal tail is then computed in subroutine 
HUPWSH. The effect of the wing lifting pressures 
acting on the nacelle area distribution is computed 
in subroutine W0NAC . 

Program WING is called by two loops in the 7,0 
level: the canard angle of attack loop and the 
pressure limiting angle of attack loop. Both 
result in changes to the local angle of attack 
distribution that cannot be handled by 
superposition. 

CALL 0VERLAY (ANLZ , 7, 2, 0) 

Configuration data in common blocks. 


See schematic, figure 2.6-2. 
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Subroutine CKLIM 


PURPOSE: Checks to see if calculated pressure coefficients 

violate limiting pressure coefficient. 

METHOD; The sum of the wing upper surface lifting pressure 

coefficient (-.5* calculated lifting pressure 
coefficient) plus the wing thickness pressure plus 
the fuselage pressure acting on the wing upper 
surface is checked against the limiting pressure 
coefficient. If the summed value, CPCHK, is more 
negative than the limit, the computed value is 
res et . 

USE; CALL CKLIM (YFR, XPC, CP, CPLIMT, PTE ST, CPT) 

INPUT: Pressure coefficient arrays in common blocks, plus: 

YFR semi-span wing station 

XPC percent chord 

CP calculated lifting pressure coefficient 

CPLIMT limiting upper surface pressure 

coefficient 

PTEST indicator variable showing limiting 

application 

CPT thickness pressure coefficient 

SUBROUTINES 

CALLED: DUBINT 


Subroutine CPTAB 

PURPOSE: To interpolate wing pressure arrays for output 

METHOD: The wing lifting pressure distributions for the 

grid system are computed in WING. CPTAB is used to 
interpolate the grid pressure distributions for 
each semi-span row at selected per cent chord 
values. The interpolated pressure coefficients are 
stored in arrays VAR (basic angle of attack) and 
VAR 1 (flat plate solution per degree alpha) . 

USE: Call CPTAB 

INPUT: Configuration data and pressure coefficients in 

common blocks. 

SUBRCUTINES 

CALLED: GGRID, TBLU1 
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Subroutine GETNCP 


PURPOSE: 

METHOD: 

use: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


To interpolate a nacelle pressure signature at a 
given wing semi-span station, ^or use in 
calculating nacelle-on-wing interference drag. 

Same as described for subroutine GETNCP for the 
near- field wave drag program. 

CALL GETNCP (YB2) 

YB2 wing semi-span fraction 


NONE 


Subroutine GETUP 

To perform double interpolation for fuselage upwash 
angle 

GETUP double-interpolates linearly for fuselage 
upwash angle in arrays TUPWC (upwash at basic 
fuselage incidence) and TUPWF (upwash per degree 
fuselage alpha) . 

CALL GETUP (XPC, YFR, UPC, UPF) 

XPC per cent chord 

YFR semi-span fraction 

UPC, UPF upwash angles, radians 


None 


Subroutine HUPWSH 
To compute downwash at tail 

HUPWSH is used to compute the local downwash acting 
along the fuselage centerline and in the plane of 
the horizontal tail. The fuselage asymmetric 
lifting pressure distribution (if any) is included 
in the wing lifting pressure definition. 

The computed downwash is stored in arrays TWT 
(basic angle of attack) and TWTE (per degree) for 
the horizontal tail. Downwash from the canard (if 
any) is contained in arrays CT and CTF and is added 
in loop D0 149. Fuselage downwash is stored in 
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USE: 

INPUT: 


SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 


arrays BX {fuselage length fraction) , BCC (downwash 
angle for basic angle of attack) and BCF (per 
degree) . 

The wing downwash can be shifted laterally to allow 
for fuselage closure effects, controlled by inputs 
FWSH and DYWH. 

CALL HUPWSH 

Configuration geometry and pressure distributions 
in common blocks 

FWSH wing downwash shift indication (0. = 

shift, 1-0 = no shift) 

DYWH downwash shift distance (y) (0- = use 

basic geometry side-of-f uselage 

increment) 


GGRID, GETR, DUBINT, TBLU1 


Subroutine P91611 

To calculate local wing slopes for grid system 

Subroutine P91611 interpolates the wing camber 
surface definition for the streamwise slopes of the 
wing grid system- Interpolation is linear spanwise 
along constant per cent chord lines, followed by 
quadratic chordwise. The resultant slopes are 
stored in array TD55DX. 

If wing twist tables or trailing edge flaps are 
input, the slopes are incremented by the 
appropriate slopes- Also, the grid element array 
containing the wing-fuselage intersection (INTN) is 
identified together with the corresponding 
fractional element (TNFR) . 

If the configuration angle of attack is not zero, 
as may be the case with limiting pressure 
calculations, all slopes are ircremented by alpha. 

If input WHUP = 1-0, the camber surface slopes are 
all zeroed. (This feature is used to generate the 
wing loading due to fuselage upwash only) . 

CALL P91611 
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INPUT: 


Configuration geometry contained in common blocks 


SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

AREA 9 
SM0G0 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


GGRID, TBLU1 


Subroutine SM00TH 

To smooth wing pressure coefficients and sum wing 
area 

SM00TH is called by program WING after all wing 
pressure coefficients have been calculated, to 
remove irregularities in the calculated values* 
Either a 9 point or 3 point smoothing equation is 
applied, depending upon the wing pressure 

calculation technique used (discussed in the theory 
document, part 1) . The wing area is also 
calculated, by summing the areas of the individual 
elements. 

CALL SM00TH (AREA9 , SM0G0) 

Configuration data contained in common blocks 
/SM0TH/, /INDEX/, and pressure data in /PC0EF) . 


(wing area of right hand wing) \/M -T . 

Smoothing technique code 

0.= 9 term smoothing 
1 - = 3 term smoothing 


None 


Subroutine W0NAC 

To calculate drag of wing lifting pressures acting 
on nacelle area distribution 

W0NAC is used to compute the thrust or drag force 
due to the wing lifting pressures acting on the 
nacelle cross-sectional distribution- For this 
calculation, the lifting pressure is broken into 
upper and lower surface halves and the proper half 
used depending upon whether the nacelles are above 
or below the wing. The pressures are transferred 
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aft along Mach lines from the wing to elemental 
frustrums describing the nacelle shape. 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 


The drag increments for both the wing pressure 
distribution at basic incidence and per degree 
alpha are computed. 

CALL W0NAC 

Configuration geometry and wing pressure field 
contained in common blocks. 

GGRID, TBLU1 

Program PUT0UT 
To print input data 

PUT0UT is used to write the input and pertinent 
program data onto the output file. 

CALL 0VERLAY (ANLZ, 7, 3, 0) 

Configuration geometry in common blocks 


None 


Program KANARD 

To calculate canard lifting pressures. 

Program KANARD calculates the canard lifting 
pressure distribution according to the same logic 
used for the wing lifting pressures. 

Call 0VERLAY (ANLZ , 7, 10, 0) 

DZDXC canard angle of attack, radians 

Configuration data in common blocks. 


GRID, GETR, DUBINT, CANPRES , CONW 
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Subroutine CAN PRES 


PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 


To sum pressure distributions over canard for lift, 
drag, and pitching moment. 

The canard pressure distributions are calculated in 
program WING, and stored in arrays TCCP (at input 
canard alpha) and TCCPF (per degree alpha) .. 
CANPRES integrates these to get lift, drag, and 
pitching moment. Drag is computed by applying the 
lifting pressure to the exposed canard slopes or 
the fuselage slopes, as appropriate. 

CALL CANPRES 

Configuration data in common blccks. 


GGRID, TBLU1 


Subroutine C0NW 


To calculate downwash from canard on wing, 
fuselage, and tail. 


Subroutine C0NW sums the downwash of the canard on 
fuselage, wing, or tail, using the same influence 
logic as the basic canard or wing solution, the 
resulting downwash arrays are: 



Influence at 

Influence per 


basic angle of 

degree alpha 


attack 


Fuselage 

FCC 

FCCF 

Wing 

TCW 

TCWF 

Hor. Tail 

TCT 

TCTF 

The downwash 

field can be 

shifted laterally 


follow the side of the fuselage, controlled by 
program input. 


CALL C0NW 


TCCP canard pressure array (basic angle 

of attack) 

TCCPF canard pressure array per degree 

alpha. 

FCSH code to shift canard downwash 

(0. = shift, -1. = no shift) 

CYCW downwash shift magnitude at wing 
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DYCH 


SUBROUTINES 

CALLED: 

PURPOSE; 

METHOD: 


(0 = use basic geometry distance 
from canard root to wing root) 
downwash shift magnitude at 
tail (0 = use basic geometry 
distance from canard root to tail 
root) . 


GGRID, GETR, TBLU1 


Program NACPF 

To calculate pressure fields acting on the wing due 
to nacelles and asymmetric fuselage volume. 

This overlay calculates the pressure fields due to 
nacelles and asymmetric fuselage volume acting on 
the wing. A schematic of the overlay is given in 
figure 2.6-3. 

The program initially calculates the asymmetric 
fuselage pressure distribution by calling 
subroutines SEGRT and SPLIT, using the fuselage 
representation described in the theory document. 

The nacelle pressure field is next calculated. A 
series of semi-span Y stations are selected, and 
the composite pressure signature due to all 
nacelles is computed; first for all nacelles below 
the wing, then for all nacelles above the wing. If 
there are nacelles both above and below the wing, 
a single nacelle lifting pressure definition is 
calculated in subroutine C0MBINE. 

The basic program format is the same as the NACPF 
subroutine described in the near- field wave drag 
program, except that the nacelle and fuselage 
interference terms are not computed in the analysis 
program version. In addition, there is an optional 
feature in the analysis program version to permit 
calculation of the nacelle pressure field at a Mach 
number other than free stream (to account for local 
Mach number effects, using input TMI^JC) . The 
resulting pressure field is afterwards referenced 
to free stream dynamic pressure- 

The following subroutines associated with NACPF are 
the same as those previously described in the near- 
field wave drag program: C0MBINE, CUT0FF, GLANCE, 
INTPLT, MERGE, NSETUP, 0RDER, PSIG. (Slightly 
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• Fuselage 
buoyancy 
field 


• Nacelle 
pressu re 
field 



Return 


Fuselage area representations 
Fuselage bouyancy pressures 

F(y) function 
Pressure signatures 

V stations for nacelle pressure field 
Nacelle setup 
F(y) for nacelles 

Pressure signatures 
Glance option 

Composite pressure 
field 


FIGURE 2.6-3.-PROGRAM NACPF SCHEMA TIC 





different name for near field program to avoid 
duplicate subroutine names) • 

USE: CALL 0VERLAY (ANLZ, 7, 4, 0) 

INPUT: Configuration geometry, Mach number contained 

in common blocks. 

SUBROUTINES 

CALLED: See schematic, figure 2.6-3. 


subroutine B0DY 

PURPOSE: To calculate Whitham F (y) function for body of 

revolution 

METHCD: Subroutine B0DY is the same as the previously 

described NB0DY subroutine in the near-field wave 
drag module, with two additional provisions: 

1) It will compute the F(y) function using the 
"smooth body" form of the F(y) equation discussed 
in reference 4, in addition to the Stieltjes 
integral equation. This provision is controlled by 
input FYLB in the calling statement. (FYLB less 
than zero uses smooth body equation) . 

2) Linear interpolation may be used (instead of 
quadratic) in fairing the body radius distribution. 
This provision is selected if FYLB is less than 
zero or greater than 9.0. 

USE: CALL B0DY (TXB, TRB, TNXB, TNCUT, TTAU, TFTAU, 

JSTO, FYLB, SWET) 

INPUT: TXB, Input body X stations and radius 

TRB, values, TNXB of each. 


TNCUT Number of body intervals in F(y) 

function 

TTAU, Y and F(y) function calculated for 

TFTAU, body, JST0 values of each 

JST0 

FYLB Calculation method code (see above) 


SUBROUTINES 

USED: TBLU1 


90 



Subr outline SEGRT 


PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


To set up area representations for below-wing and 
above wing fuselage areas 

The logic of the below-wing and above-wing fuselage 
area representations is described in the theory 
document (Part 1) . Subroutine SEGRT computes these 
area distributions using the wing- body intersection 
definition found in INTSEC. The fuselage is 
considered circular in the calculation of the two 
area distributions, resulting in: 

TRABV above-wing radius distribution 

TRBL0W below-wing radius distribution 

Alternatively, SEGRT may use input values 
of fuselage areas (if SYMM = 2.0) 

CALL SEGRT 

Configuration geometry and wing- fuselage 

intersection definition contained in common blocks 


None 


Subroutine SPLIT 

To calculate asymmetric fuselage volume pressure 
field 

Subroutine SPLIT is used to calculate the fuselage 
pressure field acting on the wing, according to 
input SYMM. If SYMM=0. , the configuration is 
considered to be mid -wing, and the input fuselage 
definition is used to calculate a symmetric (non- 
lifting) fuselage thickness pressure field. If 
SYMM> 1.0, the above-wing and below-wing fuselage 
area representations obtained in SEGRT are used to 
calculate the respective pressure fields acting on 
the wing. 

Pressure signatures due to the fuselage are 
calculated at the same X and Y stations used for 
the fuselage upwash field. The resultant fuselage 
volume pressure field is stored in common 
block/CPBASM/: 

RXUPW per cent chord array 

RYUPW semi -span percentages 
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USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 
METHOD : 

USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


FAB0VE pressure coefficients above wing 

PBEL0W pressure coefficients below wing 

CALL SPLIT 

Configuration geometry contained in common 
blocks. 

SYMM calculation code 

B0DY, PSIG, TBLU1 


Subroutine YSETUP 

To set up Y array for nacelle pressure field 
def ini t ion 

YSETUP sets up a series of Y stations located at 
each 5 per cent semi-span, plus extra stations 
located immediately inboard and outboard of each 
nacelle centerline. subroutine 0RDER is used to 
store the array (TYP) in monotcnically increasing 
fashion. 

CALL YSETUP 

Configuration geometry contained in common 
blocks. 


0RDER 


Program FINISH 

To compute fuselage forces in presence of wing 
downwash field, add in contribution of horizontal 
tail, and write out complete configuration force 
coefficient summaries. 

The wing/canard lifting pressure distribution and 
force coefficients are computed in program WING. 
Summary data from WING are passed to FINISH by 
common blocks, where the fuselage contribution in 
the wing/canard downwash field is calculated 
(subroutine FUSCF) , and the direct effects of the 
canard are added in. All coefficients are based on 
input reference geometry. 
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USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 


USE: 


A calculation loop (D0 140) then calculates the 
contribution of the horizontal tail at various 
input incidences (in subroutine HTPART) , and adds 
it to the wing-fuselage-canard data. 

A summary of the configuration force coefficients, 
nacelle-on and nacelle-off, is then computed and 
printed. This includes lift (XCL) , drag (ST0T and 
SCDN) and pitching moment (CMA and CMB) 
coefficients. Corresponding coefficients for the 
flat wing configuration are printed for reference - 

The streamwise lift distribution is then summed and 
printed (subroutine STRMWS5) . If wing lifting 
pressure coefficients at specified lift 
coefficients were requested, these are then 
calculated and printed. Finally, the wing-canard 
spanwise lift distribution is printed. 

CALL OVERLAY (ANLZ , 7, 5, 0) 

Configuration data in common blocks 


FUSCF, SUMIE, HTPART, STRMWE 


Subroutine HTPART 

To calculate the horizontal tail lifting pressure 
distribution and force coefficients 

Subroutine HTPART is used to compute lifting 
pressure distributions on the tail in the presence 
of wing-canard downwash. The equations used are 
the same as those employed for the wing; the tail 
is broken into exposed and fuselage-carry-over 
portions, and fuselage upwash and wing downwash 
added to tail alpha for the purposes of computing 
lift coefficients. 

The lifting pressure distributions are summed over 
the tail planform to get tail force coefficients, 
which are then added to the force coefficients of 
the rest of the configuration. These are passed 
back to FINISH in common blocks. 

CALL HTPART (NH, REFAR) 
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INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

use: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 


SUBROUTINES 

CALLED: 


NH Tail incidence loop index 

REFAR Reference area in program units 


GGRID, GETR, DUBINT, TBLU1 


Subroutine STRMWX 

To sum and print streamwise lift distribution 

In STRMWE, the streamwise lift distributions due to 
wing/canard, nacelles, fuselage, and horizontal 
tail are summed together and printed. The 

presentation is in fraction of total lift 
coefficient, so that the final number printed for 
the complete configuration is 1.0. 

CALL STRMW8 (REFAR, SCLN9) 


REFAR Reference area in program units 
SCLN9 Total lift coefficient 


TBLU1 


Subroutine SUMIB 

To summarize and print configuration force 
coefficients 

The configuration force coefficients, including the 
interference drag terms, are summarized for program 
FINISH in subroutine SUMIS5 . 

CALL SUMI2 (N) 

N Index variable to identify printout 

series 

1= WING 

2 — WING + FUSELAGE 

3= WING + FUSELAGE + CANARD 


None 
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Subroutine 


FUSCF 


PURPOSE: 

METHOD: 


USE: 

INPUT: 

SUBROUTINES 

CALLED: 

PURPOSE: 

METHOD: 

USE: 

INPUT: 


To calculate force coefficients of fuselage in 
presence of wing downwash field 

FUSCF repeats the slender body fuselage lift 
calculations of program FUSLGE, adding the wing 
downwash field computed in HUPWSH to the basic 
fuselage angle of attack- The fuselage is broken 
into segments and the equations given in the theory 
document (Part 1) are used to calculate the 
fuselage lift. Fuselage drag is computed by 
applying the fuselage lift distribution to the 
local mean-line slopes. 

Lift for the fuselage at basic incidence and the 
corresponding incremental (flat) fuselage at one 
degree angle of attack, plus interference drag 
terms, are all calculated. The results are stored 
in common blocks /B0DS0L/ and /EFND/: 


TBDLC 

TBDLF 

CLCND, CLFND 

CDCND, CDFND 
CMCND, CMFND 


basic incidence streamwise lift 

flat body streamwise lift (1° alpha) 

lift coefficients of basic and 

flat solutions 

drag coefficients 

pitching moment coefficients 


Data for the isolated fuselage solution (no 
wing downwash) are also printed for reference. 


CALL FUSCF 


Configuration geometry contained in common 
blocks. Isolated fuselage force data 
in /B0DS0L/. 


TBLU1 

Program L0AD9 

To read input data for lift analysis module 

Program U&AD9 reads input data (interface tape 
written by analysis subprogram of geometry module) - 

CALL OVERLAY (ANLZ, 7, 8, 0) 

See user's manual 
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SUBROUTINES 

CALLED: 


PURPOSE: 


METHOD: 


None 


Program TRNSFM9 


To convert input data to prograir units and set up 
grid system 


The wing, canard, and horizontal tail planforms are 
represented in the analysis program as a set of 
rectilinear elements, as described in the theory 
document (Part 1) . Given the number of semi -span 
element rows (FN0N) used to define the right-hand 
wing, TRNSFM interpolates the wing planform 
definition at the mid point of each semi-span row 
for leading edge and trailing edge X- value- These 
X values are converted to program scale (by means 
of RAT 10= y wing tip * BETA/FN0N) and stored "in 
arrays TXLE and TXTE, with XLE0 and XTE0 defining 
the wing centerline grid points. 


Storage of wing pressure coefficients and surface 
slopes uses a space-conserving technique: The 
leading edge of a s panwise row is stored 
immediately after the trailing edge of the adjacent 
inboard row- The index array for the row storage 
is J2W, which gives the number of elements stored 
ahead of the leading edge of a selected row- The 
same storage arrangement is used for canard and 
horizontal tail, based on indices JSC and J2H- 


The factor RATI0 is also used to scale the canard 
and horizontal tail arrays, as follows: 

Canard Tail 

Leading edge array TCXLE THXLE 

Trailing edge array TCXTE THXTE 

Centerline leading edge CXLE0 HXTE0 

Centerline trailing edge CXTE0 HXTE0 


The maximum size of the pressure coefficient arrays 
for the right hand wing is 2500 for the wing, 200 
for the canard, and 500 for the horizontal tail. 
Also, the maximum X dinension of the configuration, 
to the. most aft point on the wing, is 205 (program 
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USE: 

INPUT: 

SUBROUTINES 

CALLED: 


units) . If any of these dimensions are exceeded, 
TRNSFM rescales the program units by reducing FN0N. 


TRNSFM is also used to identify the grid elements 
associated with trailing edge flaps. Since the 
flap edges will usually not coincide with element 
edges, an approximate element array is used to 
represent the flaps and the input flap deflection 
angles are altered such that the product of flap 
area times deflection is the same for the input and 
the approximate program definition. The input flap 
area is computed in subroutine PHLAP. 


CALL 0VERLAY (ANLZ, 7, 9, 0) 


Configuration geometry read in program L0AD9 


PHLAP, GGRID, TBLU1 
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APPENDIX A 


INTERACTIVE GRAPHICS SUBROUTINES 


The graphics subroutines in the design and analysis program are 
described in this appendix, except for the LRC standard CRT 
software routines. 

Three general purpose subroutines are described first, followed by 
the subroutines associated with the individual modules; 


PURPOSE: 

METHOD; 


USE; 

INPUT; 


0UTPUT: 


Subroutine C SCALE 


Computes a plot origin and scale factor given an 
array of values. 

If the values in the input array are not equal, LRC 
subroutine ASCALE is called to compute an origin 
and scale factor. If the values in the input array 
are equal, the origin is set to that value -.5, and 
the scale factor is set to 2.0/length over which 
the data is plotted. 

CALL C SCALE (ARRAY, S,N,K,DV) 


ARRAY = 
S = 
N = 
K = 
DV = 


Array of data to be scaled. 

Length over which data will be plotted. 
Number of values in ARRAY. 

Interleave factor. (1=all points) . 
Number of divisions per inch of paper. 


ARRAY (N*K+1) = Plot origin 
ARRAY (N+K+1+K) = Plot scale factor 


SUBROUTINES 

CALLED: ASCALE 


Subroutine SHOW 
Entry SHOWI 

PURPOSE; To display an array of floating point values 

on the CDC250. 


METHCD: LRC subroutines N0TATE and NUMBER are called 

to display a variable name, equal sign and a 
number of values on the same line. When the 
display line is complete, the. vertical coord- 
inate of the line is decremented by a preset 
value. Alternate entry point SHOW I is used 
for integer values. 
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use: 

INPUT: 


OUTPUT: 

SUBR0UTINES 

CALLED: 

EXAMPLE: 

PURPOSE: 

METHCD: 

USE: 

INPUT: 


CALL SHOW (A,N,BCD,NC) 

CALL SHOWI (L,N,BCD,NC) 

A Floating point variable or array 

L Integer variable or array 

N Number of values to display 

BCD Hollerith label preceding the 

first value 

NC Number of characters in BCD 

COMMON/ PBLOK/ 

YORG Y coordinate of display line 

VORG X coordinate of first value to dis- 

play- If values to be displayed 
are negative, the negative sign 
will be positioned at VORG-CS, 
VORG+VDEL-CS and VORG+VDEL+ VDEL-CS 
YDEL Vertical distance between lines 

VDEL Horizontal distance between values 

CH Character height 

CS Spacing between characters (6/7*CS) 

THETA Angle of display for label & values 

ND Number of decimal places to display 

COMMON/PBLOK/ 

X0RG = Coordinate of first character BCD 


NOTATE, NUMBER 


R( 1 ) = 1-25 -3.79 5.86 

R ( 4 ) = 6.73 9.50 -3.17 


Subroutine SHOW3 

To display a label and three floating point 
variables on the CDC 2 50. 

LRC subroutines NOTATE and NUMBER are called 
to display a hollerith label and three values 
on the same display line. When the line has 
been displayed, the vertical coordinate of 
the line is decremented by a preset value. 

CALL SHOW3 (VI, V2, V3, BCD) 


VI = First variable to display 
V2 = Second variable to display 
V3 = Third variable to display 
BCD = Hollerith label (10 character maximum) 
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COMM0N/PBLO K / 


XORG 

YORG 

VORG 

VDEL 


YDEL 

CH 

CS 

THETA 

ND 


= X coordinate of label 

= Y coordinate of display line 

= X coordinate of VI 

= X distance between variables 

If negative, the variables are 
displayed at VORG-CS, VORG+VDEL-CS 
and VORG+VDEL+VDEL-CS 
= Y distance between display lines, 

= Character height 

= Spacing between characters (6/7+CH) 

= Angle of display for label and variables 
= Number of decimal places to display. 


SUBROUTINES 

CALLED: NOTATE, NUMBER 

EXAMPLE: 

NAC1 -57.254 9.000 5.790 

NAC2 69.874 16.785 -9.103 


GEOMETRY MODULE 


Configuration geometry is displayed and/or edited by program 
DISGE0M, a secondary overlay in the geometry module. 


PURPOSE: 

METHOD: 


USE: 


Program DISGE0M 

To display and/or edit configuration geometry 

A description of the display capability is 
presented in Appendix A of the user's manual. 

CALL OVERLAY (GEOM, 5, 5, 0) 


Subroutine ALTER 

PURPOSE: To alter the wing camber surface shape to match a 

new trailing edge definition. 

METHOD: The camber definition at each input airfoil is 

rotated about the airfoil leading edge point until 
the trailing edge point coincides with the new 
trailing edge definition. The new camber surface 
definition is stored in temporary array B80RD 
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USE: 


INPUT: 


OUTPUT: 

PURPOSE: 

METHOD: 


USE: 

INPUT: 

OUTPUT: 


CALL ALTER {2TE) 

ZTE = Array of new trailing edge values. 
COMMON/TEMP/ 

DECODE = +, 2TE array consists of camber values 
DECODE = 2TE array consists of camber values 
+ the 2 value of the leading edge 


COMMON /WING/ 

T20RD = Original camber definition 


COMMON/TEMP/ 

B20ED = Altered camber definition 


Subroutine PLTSI2 

To compute the configuration minimum and maximum X 
and Y coordinates, and the fuselage minimum and 
maximum 2 coordinates. A plot scale factor is also 
output . 

Each configuration component is analyzed as to its 
coordinate values. The minimum and maximum values 
are stored in common. If the range of data in the 
X direction is greater than that in Y, SCALE is 
computed as: 

X MAX -XMIN 

10.0 

If the range of data in the Y direction is greater 
than that in X, SCALE is computed as: 

YMAX-Y MIN 

7.0 


CALL PLTSI2 (SCALE) 

All configuration geometry in COMMON. 

SCALE = Plot scale factor 
COMMON /0VL1 / 

XMIN = Minimum X value of configuration 

XMAX = Maximum X value of configuration 

YMIN = Minimum Y value of configuration 

YMAX = Maximum Y value of configuration 

2MIN = Minimum « value of fuselage 

2MAX = Maximum 2 value of fuselage 
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FAR-FIELD WAVE DRAG MODULE 


The graphics subroutine (DIS080) in the far- field wave drag module 
is located in program 0UT, overlay (FFWD,3,6) . 

Subroutine 0UT 


PURPOSE: 

METHOD: 


USE 

INPUT: 


To display area plots and drag summary 

The CRT is used to display far- fie Id wave drag 
module results as described in the user's manual. 
Appendix A. 

CALL DIS080 (S # B, B0, C, RC, N) 

S Array of X values 

B Array of fuselage areas 

B0 Array of optimum fuselage areas 

C Array of overall configuration areas 

RC Array of restrained areas 

N Number of values in input arrays 


NEAR -FI ELD WAVE DRAG MODULE 


The graphics program (DISP916) is called from the primary level of 
the near-field wave drag overlay. 


PURPOSE: 

METHOD: 

USE: 
INPUT : 


Program DISP916 

To display near-field wave drag module results 

DISP916 displays summary results described in 
the user's manual. Appendix A 

CALL 0VERLAY (NFWD, 6, 3, 0} 

Summary data in common blocks 

SKIN FRICTION DRAG MODULE 


The graphics subroutine (DISP158) in the skin friction module is 
called from subroutine DRAG. 
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Subroutine DISP158 


PURPOSE: To display skin friction module results 

METHOD: DISP158 displays summary results described in the 

user’s manual. Appendix A. 


USE: CALL DISP158 

INPUT: Summary data in common blocks 


WING DESIGN MODULE 


Display and editing in the wing design module is performed by 
three graphics programs called from the primary overlay. 


PURPOSE: 

METHOD: 


USE: 

INPUT: 


Program BUCKETP 

To display drag-due -to-lift bucket plot 

BUCKETP is used to display the K^. versus Cmo plot 
plus design point solutions described in Appendix 
A of the user’s manual. 

CALL 0VERLAY (WDEB, 10 , 4, 0) 

Data in common block/ BUCKET/ 


PURPOSE: 

METHOD: 


USE: 


Program EDBUCK 

To permit editing of design point variables 

EDBUCK allows user to edit wing design program 
variables Cmo, CLDBIN, RESTART, and C0NSTR(1) 
through C0NSTR (4) . In addition it allows the user 
to execute the next design case or calculate the 
edited design point, as described in Appendix A of 
the user’s manual. 

CALL 0VERLAY (WDEB , 10, 5 , 0 ) 


Program ST0P0P 

PURPOSE: To permit termination of wing design program 

cases 

METHOD: ST0P0P allows user to exit from a series of wing 
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design cases, as described in Appendix A of "the 
user's manual. 

USE: CALL 0VERLAY (WDE2 , 1 0, 6, 0) 

LIFT ANALYSIS MODULE 


The graphics displays in the analysis 
overlays called from the primary overlay, 
called from program FINISH. 


module consist of two 
plus two subroutines 


PURPOSE: 

METHOD: 


USE: 

INPUT: 


Program DISTWST 

To display wing twist and permit editing of twist 
and several execution codes. 

DISTWST displays and permits editing of wing twist 
array, plus canard angles of attack, SYMM, WHUP, 
and ANYBOD. The display presentation is described 
in Appendix A of the user 1 s manual. 

CALL 0VERLAY (ANL2, 7, 6, 0) 

Input variables in common blocks 


Program DISUPWS 


PURPOSE: 

METHOD: 

USE: 

INPUT: 


To display fuselage upwash and wing pressure 
coefficient data 

DISUPWS provides display options for calculated 
fuselage upwash or wing pressure coefficient data, 
as described in Appendix A of the user's manual. 

CALL 0VERLAY (ANL£, 7, 7, 0) 

Data in common blocks /UPWSH/and/CPBUPWI/. 


Subroutine D I STAB 

PURPOSE: To display analysis module force coefficient 

summary and permit editing of horizontal tail angle 
of attack. 

METHOD: DISTAB provides display of analysis program results 

as described in Appendix A of the user's manual. 
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USE: 


CALL DISTAB (NH, HTALP) 


INPUT: 


NH 

HTALP 


Horizontal tail angle of attack 
loop index 
Edited tail alpha 


PURPOSE: 


METHOD: 


USE: 

INPUT: 


Subroutine EXL00P 

To display and permit editing of variables 
within DO loops. 

EXL00P displays current values of Mach number, 
configuration alpha, and canard alpha. It 
then permits editing of the next value to 
change within the cycle, as described in 
Appendix A of the User’s Manual. 

CALL EXL00P 

Configuration data in common blocks 
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